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A STUDY  ON  THE  PERFORMANCE  OF  RESIDENTIAL  BOILERS  FOR  SPACE  AND 

DOMESTIC  HOT  HATER  HEATING 

ABSTRACT 

A residential  boiler  for  space  heating  and  domestic  hot  water  heating  was 
studied  by  conducting  laboratory  tests  and  computer  simulations.  A clam- 
shell, wet-base,  oil-fired,  residential  boiler  with  a tankless  coil  for 
heating  domestic  water  was  selected  for  this  research  project. 

The  purpose  of  this  study  was  to  develop  a method  for  evaluating  the 
performance  of  an  integrated  space  and  water  heating  appliance.  Based  upon 
laboratory  tests,  a computer  model  was  developed  and  used  with  the  HVACSIM^ 
building  system  simulation  program  to  simulate  the  operation  of  the  integrated 
appliance . 

The  model  was  verified  for  heat-up,  cool-down,  cyclic,  and  standby  modes  of 
operation,  along  with  various  domestic  hot  water  draw  cycles.  Using  the 
verified  model,  computer  simulations  were  carried  out  for  both  summer  and 
winter  operations  of  the  appliance.  As  a result  of  these  simulation 
studies,  a simple  method  for  determining  the  combined,  seasonal  efficiency 
of  Type  I appliance,  whose  primary  design  function  is  space  heating  and 
secondary  function  is  domestic  water  heating,  is  presented. 
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1.  INTRODUCTION 


Although  boilers  with  tankless  coils  have  been  around  for  a long  time,  other 
types  of  residential  integrated  appliances  designed  for  both  space  and 
domestic  hot  water  heating  have  recently  emerged  in  the  market  place. 

Because  of  the  newness  of  these  combined  appliances,  a new  method  for  rating 
their  performance  has  become  necessary.  Responding  to  this  need,  the  ASHRAE 
SPC-124P  committee  has  developed  a draft  standard  [1]  entitled  "A  Method  of 
Testing  for  Rating  Combination  Space  Heating/Water  Heating  Appliances"  and 
submitted  it  for  public  review. 

Previously,  Subherwal  [2],  Pietsch  [3],  and  Nordstrom  and  Fuller  [4]  discussed 
the  performance  rating  of  combination  space  heating/water  heating  appliances 
in  several  papers.  In  an  effort  to  review  the  ASHRAE' s proposed  test  method, 
laboratory  tests  and  computer  simulations  were  performed  at  the  National 
Institute  of  Standards  and  Technology  (formerly  National  Bureau  of  Standards) 
on  a clam-shell,  wet-base,  oil-fired,  residential  boiler  with  a tankless 
domestic  water  heating  coil.  Based  upon  the  commercial  steam  boiler  model 
by  Chi  [5],  a residential  boiler  computer  model  was  developed  and  is 
discussed  in  this  report.  Some  of  the  information  from  a commercial  boiler 
model,  which  was  developed  by  Chi,  Chern,  and  Didion  [6],  was  also  incorporated 
in  this  boiler  model. 

The  computer  model  was  then  used  to  simulate  the  operation  of  the  selected 
appliance  by  incorporating  it  in  the  HVACSIM^  program,  which  is  a dynamic 
building  system  simulation  program  created  at  NIST  [7,8].  The  computer 
model  was  verified  with  laboratory  test  data  for  heat-up,  cool-down,  cyclic. 
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and  standby  modes  of  operation,  along  with  various  domestic  hot  water  draw 
cycles . 

A series  of  computer  simulations  were  performed  for  domestic  water  heating 
only  (summer  operation)  and  for  combined  space  and  domestic  water  heating 
(winter  operation)  . A family  of  combined  efficiency  cutrves  was  obtained 
as  a function  of  the  space  load  factor  for  different  domestic  water  loads. 

An  analysis  of  the  behavior  of  the  selected  boiler  with  a tankless  coil 
leads  to  a new  method  for  determining  the  combined,  seasonal  efficiency  of 
appliances  whose  primary  and  secondary  functions  are  space  and  domestic 
water  heating,  respectively. 

2.  LABORATORY  TESTINGS 

2.1.  Test  Setup 

The  integrated  appliance  under  study  was  a residential  wet-base  hot  water 
boiler  with  five  sectional  clam-shell,  cast-iron  heat  exchangers.  One  of 
the  heat  exchanger  sections  contained  a finned,  copper  coil  for  domestic  hot 
water  heating.  The  boiler's  DOE  rated  heating  capacity  was  158  kBtu/h  and 
the  firing  rate  of  the  oil  burner  was  1.35  gal/h.  The  boiler  used  a No. 2 
fuel  oil.  The  external  view  and  the  heat  exchanger  arrangement  of  the 
boiler  are  shown  in  Figures  1 and  2,  respectively. 

Figure  3 shows  the  overall  sketch  of  test  setup,  excluding  some  of  the  data 
acquisition  instruments.  A six- inch  diameter  stack  was  directly  connected 
to  the  boiler  top,  and  an  1/125  HP  circulating  water  pump  for  space  heating 
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Figure  3 Test  set-up  of  the  boiler  with  a tankless  coll 


loop  was  located  near  the  return  port  to  the  boiler.  Two  one-inch  pipes 
were  also  connected  to  the  inlet  and  the  outlet  ports  of  the  domestic  water 
heating  coil.  The  majority  of  water  piping  was  insulated  with  1/2 -inch 
insulation.  A barometric  damper  for  draft  control  was  located  in  the 
stack. 

Two  load  tanks  were  used  to  simulate  the  building  load.  Cold  water  from 
the  tap  cooled  down  the  circulating  boiler  water  in  the  main  load  tank. 

An  auxiliary  load  tank  served  to  increase  the  heat  removal  capacity  of  the 
cooling  water.  A pump  circulated  the  cooling  water  through  both  load  tanks. 

Only  the  main  load  tank  contained  a heat  exchanger  for  cooling  the  circulating 
boiler  water. 

A controller,  supplied  by  the  manufacturer,  controlled  the  operation  of  the 
oil  burner.  The  input  signal  to  the  controller  was  provided  by  the  thermocouple 
located  in  the  clam-shell  section  of  the  heat  exchanger  which  contained  the 
tankless  coil.  The  controller  also  controlled  the  water  circulating  pump 
depending  upon  the  on/off  condition  of  a thermostat  that,  in  this  case, 
sensed  the  load  tank  water  temperature.  An  override  switch  was  installed  to 
bypass  this  thermostat  and  permit  manual  control  of  the  pump. 

A number  of  thermocouples  were  installed  at  various  locations  in  the  test 
apparatus.  These  thermocouples  were  connected  to  a data  acquisition/control 
instrument  that  was  connected  to  a personal  computer.  The  test  was  automated 
using  off-the-shelf  data  acquisition  and  control  application  software.  A 
gas  analyzer  was  used  to  measure  the  CO2  concentration  in  the  flue  gas  of 
the  boiler. 
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The  water  flow  rate  for  the  space  heating  function  was  set  by  a manually 
operated  valve,  while  the  flow  rate  for  domestic  hot  water  remained  constant 
during  all  water  draws.  Two  water  flow  meters  with  electric  pulse  generators 
were  used  for  measuring  the  amount  of  water  flow  through  the  boiler  and  the 
tankless  coil.  An  in-line  solenoid  valve,  controlled  by  the  computer  based 
data  acquisition/control  system,  controlled  the  flow  of  domestic  hot  water. 

2.2.  Test  Procedures 

Laboratory  testings  were  performed  to  provide  input  data  for  the  computer 
model  and  to  verify  the  model.  Experimental  work  involved  heat-up  and 
cool-down  tests,  a steady-state  test,  and  tests  to  determine  the  effects  on 
efficiency  of  various  space  loads,  water  loads,  and  combined  space  and 
domestic  water  loads. 

Prior  to  a test,  the  data  acquisition  sequence  and  control  sequence  were 
programmed  into  the  computer.  Measured  data  was  automatically  stored  on  the 
hard  disk  in  the  personal  computer.  Measurements  of  oil  consumption  and 
electric  energy  use  were,  however,  performed  manually.  Oil  consumption  was 
determined  by  recording  the  weight  of  the  oil  container  before,  during,  and 
after  a test.  Electric  energy  consumption  was  measured  manually  reading  a 
watt-hour  meter. 

The  data  sampling  period  was  limited  by  the  size  of  data  file  that  could  be 
stored  on  the  hard  disk.  In  order  to  keep  the  data  file  manageable,  two 
periods  were  used.  A small  sampling  period  was  used  for  fast  changing 
situations,  while  a large  one  was  used  for  slow  changes. 
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During  a test,  a color  monitor  displayed  plots  of  selected  variables  with 
respect  to  time.  The  displayed  plots  revealed  very  valuable  information  on 
the  status  of  the  test.  Figure  4 is  a sample  screen  dump  to  a printer.  The 
numeric  value  appearing  inside  a rectangular  box  is  the  totalized  value  of 
water  flow  in  gallons  for  this  particular  test. 

Heat-up  and  cool-down  tests  were  made  without  any  external  water  flow,  i.e. 
no  space  heating  load  or  domestic  hot  water  load  was  imposed.  The  oil  burner 
of  the  boiler  stayed  on  until  the  boiler  water  temperature  reached  its  cutoff 
point.  Since  no  external  water  flow  was  allowed,  the  boiler  water  cooled 
down  very  slowly. 

During  steady-state  tests,  the  space  heating  water  circulated  continuously 
through  the  boiler  and  the  heat  exchanger  of  the  main  load  tank.  The  water 
flow  rate  was  adjusted  to  make  the  burner  run  continuously  without  causing 
significant  variation  of  the  boiler  water  temperature.  It  was  found,  however, 
that  it  was  very  difficult  to  make  such  an  adjustment  due  to  variations  in 
boiler  and  load  tank  temperatures.  Because  of  this,  steady- state  experiments 
usually  lasted  less  than  30  minutes.  No  domestic  hot  water  was  drawn  during 
these  tests. 

The  performance  rating  test  was  carried  out  following  the  test  procedure 
given  by  the  ASHRAE/ANSI  Standard  for  rating  the  performance  of  residential 
boilers/furnaces  [9] . The  quantities  required  by  the  Standard  were  measured 
and  then  used  to  calculate  the  boiler  seasonal  fuel  utilization  efficiency 
and  steady-state  efficiency. 
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Figure  A 


Sample  screen  dump  of  the  displayed  plots  on  a monitor 


Space  heating  loads  were  simulated  through  the  use  of  a load  tank  previously 
described.  The  water  temperature  within  the  load  tank  was  controlled  by 
means  of  the  thermostatically  controlled  water  pump.  The  thermostat's  upper 
and  lower  limits  were  set  before  each  test.  No  domestic  water  was  drawn 
during  a space  load  test  and  the  controller  supplied  with  the  boiler 
controlled  the  burner  operation.  Using  the  water  flow  rate  and  the  temperature 
difference  between  the  boiler's  supply  and  the  return  water,  the  amount  of 
energy  delivered  to  the  load  tank  was  computed.  The  effect  of  the  water 
pump  cycling  rate  on  the  boiler's  energy  consumption  was  also  investigated. 

Figure  5 illustrates  how  the  burner  and  pump  on/off  status  and  the  stack  gas 
temperature  changes  with  respect  to  time  for  a typical  space  load  simulation 
test.  Due  to  the  limitation  of  the  load  tank  cooling  capacity,  the  cycle 
rate  could  not  be  increased  over  a certain  limit.  The  duration  of  a typical 
test  with  repeatable  cycles  was  usually  between  2 and  3 hours. 

Domestic  hot  water  load  simulation  tests  were  carried  out  according  to  the 
ASHRAE  12AP  proposed  Standard  [1].  Even  though  the  measured  first  hour 
rating  was  138  gal,  the  maximum  allowable  total  daily  draw  of  120  gal 
specified  in  the  proposed  Standard  was  used  in  all  draw  tests.  After  the 
18-hour  standby  period,  six  equal  draws  of  domestic  hot  water  were  imposed 
at  the  beginning  of  each  hour.  The  sampling  rate  of  data  collection  was  one 
scan  per  15  seconds. 

During  the  water  load  simulation  tests,  no  space  load  was  applied.  The  lower 
and  upper  setpoints  of  the  boiler  controller  were  set  to  190° F (87.78°C)  and 


10 


o STACK 


e 8 

TiMc  iS) 
SUPINER 


PUMP 


Figure  5 The  burner  and  pump  on/off  status  and  the  stack  gas  temperature 
of  a typical  space  load  simulation  test 
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210°F  (98.89°C),  respectively.  The  flow  rate  of  domestic  hot  water  was  fixed 
as  assigned  by  the  Standard  at  3.0  gal/min. 

In  the  18-hour  standby  period,  the  boiler  water  temperature  was  maintained 
between  these  setpoints  by  operating  the  burner  to  compensate  the  heat 
losses  due  to  the  stack  gas  flow,  heat  flow  through  the  boiler  jacket,  and 
heat  conduction  to  pipes. 

By  combining  space  load  and  domestic  water  load,  combined  load  tests  were 
performed.  The  simulation  tests  were,  however,  restricted  to  low  space  loads, 
due  to  the  limited  cooling  capacity  of  the  load  tanks.  The  water  pump  circulated 
the  space  heating  water  with  a constant  flow  rate,  when  space  heating  was 
demanded.  To  simulate  the  domestic  water  load  simulation  tests,  one-sixth 
of  the  daily  usage  of  hot  water  was  drawn  at  a fixed  rate  at  the  beginning 
of  each  hour  of  the  six-hour  draw  period.  The  combination  of  two  space 
heating  and  domestic  water  heating  loads  resulted  in  complex  cycles  as,  for 
example,  shown  in  Figure  6.  The  status  of  the  burner  and  circulating  pump 
and  the  stack  gas  temperature  are  shown  in  this  figure. 

3.  COMPUTER  SIMULATIONS 

The  HVACSIM^  program  was  used  for  computer  simulations  of  the  boiler  with  a 
tankless  coil.  Component  models  consisting  of  a water  boiler  model,  a 
simple  heating  coil  model,  and  an  algorithm  for  boiler  control,  were 
developed  to  be  compatible  with  the  HVACSIM^  program.  These  component 
models  were  connected  to  each  other  to  form  a model  of  the  boiler  with  a 
tankless  coil. 
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Figure  6 The  status  of  the  burner  and  circulating  pump  and  the  stack  gas 
temperature  of  a combined  load  test 


The  boiler  model  was  an  empirical  model  that  required  reasonably  good 
initial  input  data  based  on  the  actual  data.  The  input  data  were  prepared  from 
laboratory  tests  and  the  boiler  configuration.  As  will  be  discussed  later 
in  this  report,  a number  of  simplifying  assumptions  were  made  in  modelling 
the  boiler/tankless  coil  system. 

3 . 1 Computer  Model 

3.1.1.  Water  Boiler  Component  Model 

Figure  7 shows  a cross-sectional  view  of  the  boiler  that  was  modelled.  As 
seen  in  the  figure,  a domestic  water  heating  coil  is  located  in  the  right- 
hand  side  of  the  heat  exchanger.  For  modelling  purpose,  the  boiler  and  the 
coil  were  considered  as  separate  component  models. 

The  empirical,  residential,  fossil  fuel-fired,  water  boiler  model  was  based 
on  the  simplified  commercial  boiler  model  by  Chi  [5].  The  schematic 
diagram  of  the  boiler  model  is  depicted  in  Figure  8.  The  supply  water 
temperature  to  the  load  was  assumed  to  be  the  same  as  the  boiler  water 
temperature.  The  boiler  water  temperature,  in  turn,  was  assumed  to  be 
uniform  inside  the  boiler.  These  assumptions  were  made  to  simplify  the 
modelling  task. 

As  shown  in  Figure  8,  the  boiler  can  be  divided  into  the  following  five 
sections:  combustion  gas  product,  fire-box  wall,  heat  exchanger  wall,  boiler 
water,  and  boiler  jacket.  The  heat  transfer  phenomena  within  the  boiler  are 
different  during  the  on  and  off  periods  of  the  burner. 


14 


Figure  7 


Stack  gas 


Jacket 


The  cross-sectional  view  of  the  boiler 


^ 9pw 


I’stk 

HXwan 

Qhx 

Boiler 

Jacket 

Gas 

■*"  water 

°i 

^lahx 

V 

Rre-box 

wall 

J 

Room  air 


'ra 


^input 


Figure  8 


A schematic  diagram  of  the  modelled  boiler  sections 


3. 1.1.1  Gas-side  Heat  Transfer 


Referring  to  Figure  8,  a heat  balance  on  the  gas-side  can  be  represented 
by: 

Qinput"Qlat”Qf"*^x"Qstk~®-  (1) 

In  this  equation,  Qinput  fuel  input  rate  which  is  computed  from  the 

mass  flow  rate  of  fuel,  , and  its  higher  heating  value  of  fuel,  HHV, 

using 


Qinput  = %u.l  HHV.  (2) 

• 

The  quantity  latent  heat  loss  calculated  using: 

Qiat  ~ 2442.0  %ual^tf*  (3) 

where  R^tf  ratio  of  mass  of  water  in  the  fuel  to  the  total  mass  of 

the  fuel  and  the  constant  2442.0  (kJ/kg)  is  latent  heat  for  evaporation  of  water. 

_ • 

The  term  is  the  heat  flow  rate  from  the  combustion  gas/draft  air  to  the  boiler 
water  through  the  wall  of  the  fire-box.  The  term  0^,3^  denotes  the  heat 
transfer  rate  between  the  gas  and  the  boiler  water  through  the  heat-exchanger. 
Once  these  four  heat  flow  rates  are  determined,  the  heat  loss  through  the 
stack  is  Q,tk  obtained  using  equation  (1). 

3. 1.1. 2 Fire-box  Wall  Heat  Transfer 

Heat  transfer  through  the  fire-box  wall  during  the  burner  on-period  was 
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calcuiated  using  the  well-stirred  combustion  chamber  theory  [10].  The 
radiative  and  convection  heat  transfer  from  combustion  gas  to  a sink 
surface  is  given  by: 


g , a V e 


A _ T 

■‘•k  s > 


(4) 


where  a is  the  Stefan-Boltzmann  constant  (5.67*10”^^).  The  absolute  temperature 
Tg  ave  average  gas  absolute  temperature  of  the  fire-box  exit  gas 

temperature  and  the  absolute  adiabatic  flame  temperature,  T^f.  The  absolute 
adiabatic  flame  temperature  is  given  by: 


(5) 


where  T 


ra»  ^p,  product  ^tf  boiler  room  air  temperature,  the  specific 


heat  of  the  combustion  products  and  the  mass  ratio  of  combustion  products  to 
fuel,  respectively.  The  absolute  temperature  of  the  sink  surface,  Tj^  3 , is 
assumed  to  be  the  same  as  the  absolute  boiler  water  temperature.  The 
quantity  is  an  effective  emissivity,  the  sum  of  radiative  and  convective 

heat  transfer  parts  and  given  by: 


A,,  = G,  (1  - K3)/(1  - K*)  + 2 A,f  h,/[a  (T, 


g , a V e 


(6) 


where 


(7) 


^ ~ ^k  s / ^g , ave  • 
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In  the  above  equation,  A^.  ^ and  A^.  ^ are  the  surface  areas  of  convective  heat 
transfer  and  radiative  heat  transfer,  respectively;  h^.  is  the  convective 


heat  transfer  coefficient:;  €gas  ^sink  emissivity  and  the 

sink  surface  emissivity  of  the  boiler  fire-box;  and  Cg  is  the  area  ratio  of 
boiler  fire-box  effective  radiation  heat  transfer  area  to  the  total  boiler 
fire-box  area. 


Using  the  Nusselt  number,  Nu,  expressed  in  terms  of  Reynolds  number.  Re,  and 
Prantl  number,  Pr,  the  convective  heat  transfer  coefficient,  h^  can  be  evaluated. 

Nu  = 0.023  Re°-®  Pr° • ^ for  Re  > 2000,  (8) 

Nu  = 3.66  for  Re  < 2000,  (9) 

and  hg  = Nu  k / dj^  , (10) 

where  k is  the  thermal  conductivity,  and  the  hydraulic  diameter,  dj^  , can  be 
represented  in  terms  of  the  volume  of  fire-box  gas,  and  the  gas  path 

length,  Lgpf,  using: 

= 1 4 / (,  L„,)  (11) 

During  the  on-period  of  the  boiler  burner,  the  mass  flow  rate  of  combustion 
product,  mg,on>  given  by: 

, on  = “fuel  ^tf 
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The  Reynolds  number  and  Prantl  number  used  in  equation  (8)  are  given  by 


R®  = “g.  on  dh  /(>  Across  ).  and  (13) 

Pr  = /i  Cp  / k,  (14) 


where  n is  the  dynamic  viscosity,  Cp  is  the  specific  heat,  and 
the  cross  sectional  area  through  which  a fluid  passes. 


During  the  off -period  of  the  burner,  the  mass  flow  rate  of  the  draft  air, 

• 

off>  ®an  be  obtained  using  an  equation  similar  to  the  one  used  in 
ASHRAE/ANSI  Standards  [9],  or  in  the  report  by  Kelly,  Chi,  and  Kuklewicz  [11]: 

- T,.)/(T„,,..  - T„)l" 

* KT.tk,..  + 273.15)/(T,,,  + 273.15)]‘  »\  (15) 

where  ,,  is  the  stack  gas  temperature  at  steady  state. 

The  effectiveness  method  for  a compact  heat  exchanger  [12,13]  is  applied  to 
calculate  the  off-period  convective  heat  transfer,  off- 

Qf.otf  = Cp.  (T,.  - Tp.)  (1  - e-»”),  (16) 

where  NTU  = A.,  h.  / (Cp,  (12) 

NTU  is  the  number  of  heat  transfer  units,  Cp^  is  the  specific  heat  of  air,  and 
Tj,^  is  the  boiler  water  temperature.  The  exit  gas  temperature  of  the  fire- 
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box  during  the  off-period,  is  obtained  using  equation  (16): 


T = T 

■^ex.g.f  -^ra 


Qf.off  / 


pa  *^8  , o f f 


) 


(18) 


During  the  on-period,  the  exit  gas  temperature  is  estimated  iteratively  using 
Newton's  method. 


3. 1.1. 3 Heat  Exchanger  Heat  Transfer 

The  number  of  heat  transfer  units  between  the  gas  and  the  boiler  water  at 

the  full  load  condition,  NTUf , can  be  determined  from  the  boiler  test  data  using 

a semi-empirical  equation  [14]: 

NTU,=  ln(l  - - Q,  - Q,)/[C„  (19) 


where  Qj  , Cp^  , steady-state  efficiency,  the 

jacket  heat  loss  rate,  the  specific  heat  of  gas,  the  gas  temperature  at  the 
inlet  of  the  heat  exchanger,  and  the  supply  water  temperature,  respectively. 
With  the  boiler  heat  exchanger  heat  transfer  number  at  full  load  given  by 
equation  (19),  the  number  of  heat  transfer  units  at  part  load,  NTU  on,  can 
be  evaluated  by  scaling  as  follows: 


NTU,on=NTU,  (m^ , ^ /m, , ^ p )° • ^ (m* . on . f /^g , on . p ^ 

*(Cp,on.f/Cp,on.p)°  ® (kon,pAon,f)°  ® . 


(20) 


where  the  subscript  p indicates  part-load  condition  and  f indicates  full 
load. 
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The  heat  flow  rate  during  the  on-period  between  the  combustion  gas  products 


and  the  boiler  water,  Q^x  on*  obtained  using  equation  (20) 


Qhx.on  = - T,,)  (1  - 


(21) 


Similarly,  the  heat  flow  rate  during  the  off-period  is 


Qh: 


off  g 


off  ™g , o f f 


(Tin. 


h X 


- T,.)  (1  - e 


■NTU .off 


). 


(22) 


The  off -period  number  of  heat  transfer  units,  NTU, off,  is  determined  based 
on  gas  properties  during  the  off-period.  When  the  gas  flow  rate  is  very 
small,  equation  (20)  becomes  infinite,  and  the  following  expression  is 
instead  used: 


^x.off  “ ^pg.off  “g.off  (T^n.hx  ■ "^sw)-  (23) 

3. 1.1. 4 Boiler  Water  Heat  Transfer 

Dynamic  changes  in  the  boiler  water  temperature,  the  most  important  quanti  ty 
in  the  boiler  model,  are  considered  using  an  ordinary  differential  equation 
as  follows : 

Cp,  M.  (dTb,/dt)  = Cp,  i,  (T„  - T^,)  + Q,.,  (24) 

For  this  relation,  the  boiler  water  temperature  is  assumed  to  be  the  same  as 
the  supply  water  temperature  to  the  space  load.  Cp^  is  the  specific  heat  of 
the  boiler  water,  Mw  is  the  mass  of  boiler  water,  m„  is  the  circulating 
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water  flow  rate,  Tj.^  is  the  return  water  temperature,  and  is  the  heat 

gain  rate  of  the  boiler  water  given  by 

Qss  = Qf  + Qhx  - Qj  - Qcw.  (25) 

where  is  the  heat  flow  rate  from  the  boiler  water  to  the  tankless  coil  water. 

Considering  the  temperature  lag  due  to  the  thermal  mass  of  the  heat  exchanger 
and  boiler  water,  the  heat  gain  at  the  current  time  can  be  replaced  by  the 
heat  gain  at  the  previous  time  step.  In  addition,  the  capacitance  of  the 
boiler  water,  Cp„  , can  be  given  as  the  nominal  time  constant  times  a 
constant,  that  is  determined  empirically.  Equation  (24)  can  thus  be 

written  as 

Ic  f w w ) ~ w C^r  w ' w ^ Qs  s , - 1 > 

where 

Q8»,-i(t)  = QssCt  - At), 

is  an  empirically  determined  integration  constant,  and 
is  the  nominal  boiler  water  time  constant. 

3. 1.1. 5 Boiler  Jacket  Heat  Transfer 

Heat  transfer  through  the  boiler  jacket  from  the  boiler  water  to  the 
ambient  air  is  calculated  from  the  overall  heat  transfer  coefficient,  Uj  , 
the  jacket  area,  Aj  , and  the  temperature  difference  across  the  jacket,  AT. 

Qj  = Aj  Uj  AT  = Aj  Uj  (Tb,  - T,.)  (28) 


(26) 

(27) 
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since  the  mass  of  the  jacket  is  much  lower  than  the  mass  of  the  boiler 
water  or  that  of  the  heat  exchanger  , instantaneous  thermal  response  is 
assumed  in  the  equation  above. 

3. 1.1. 6 Stack  Gas  Temperature 


Thermal  properties  of  the  gas  in  the  heat  exchanger  such  as  /i,  Cp  , and  k, 
are  determined  using  the  average  temperature  of  the  inlet  and  outlet  gas 
temperature  of  the  heat  exchanger.  The  inlet  gas  temperature  is  the  exit 
gas  temperature  of  the  fire-box,  but  the  outlet  gas  temperature  is  the 
stack  gas  temperature  as  shown  in  Figure  8.  Moreover,  the  off -period  mass 
flow  rate  of  the  gas  depends  upon  the  stack  gas  temperature  (see  equation  15). 

Separate  differential  equations  for  the  stack  gas  temperature  are  considered 
for  the  on-period  and  the  off-period,  since  the  patterns  of  rising  and  decay 
of  the  stack  gas  temperature  are  usually  different  in  each  case. 

During  the  on-period. 


^cfg.on  "^stk 


(29) 


and  during  the  off -period, 

^efg.off  ’'g  (T.tk/^^^)  ^Btk  ~ ^bw.taf  + (1  " w)  Tj.  ^ , (30) 

where  is  a nominal  gas  time  constant^  and  Icfg.on  ^cfg.off 

integration  constants.  Appropriate  integration  constants  and  the  real  time 
constants,  for  the  on-period  and  for  the  off-period,  can  be  determined 
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based  on  laboratory  test  results.  The  quantity  33  is  the  steady-state 

stack  gas  temperature,  Tj,^  ref  a reference  boiler  water  temperature 
corresponding  to  the  stat:k  gas  temperature  at  an  equilibrium  state,  and  w is 
a weighting  factor  that  is  used  to  obtain  a good  fit  to  the  measured  stack 
gas  temperature  decay  curve.  The  values  of  Tj,^  ref  used  in  most  of 

computer  simulations  were  97. 5° C and  0.8  respectively.  These  assigned 
values  can  be  changed  depending  upon  the  characteristics  of 
a boiler  of  interest. 

It  should  be  noted  that  the  stack  gas  temperature  specified  by  equations 
(29)  and  (30)  bound  the  heat  transfer  rates  of  the  heat  exchanger  during 
the  on-  and  off -periods. 

3.1.2.  Domestic  Water  Heating  Coil  Model 

The  coil  model  for  domestic  water  heating  was  developed.  This  model  was  an 
extension  of  a simplified  approach  for  calculating  the  heat  transfer 
resulting  from  flow  through  a pipe  with  a constant  surface  temperature, 

^surf  c model,  the  boiler  water  temperature  is  assumed  to  be 

uniform  everywhere  inside  the  boiler,  and  the  capacitance  of  the  coil  wall 
is  neglected.  Properties  of  the  water  in  the  coil  are  evaluated  at  the 
average  temperature  of  the  coil  inlet  and  outlet  temperatures. 

The  convective  heat  transfer  coefficient  between  the  coil  wall  and  the  coil 
water,  h^.^,  is  calculated  in  a similar  manner  as  given  by  equation  (10). 

The  value  of  h^  is  obtained  for  given  values  of  the  pipe  diameter,  d^ , the 
pipe  length,  Ig  , the  inlet  temperature,  g^,  and  the  mass  flow  rate  of 
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coil  water,  The  outlet  temperature  of  the  coil  water,  cw> 

the  heat  transfer  rate,  can  be  computed  from: 


T = T - T 

out.cw  Burf.c  ' surf.c 


(31) 


Qc  w w u r f , c m > 


(32) 


where 


~ (^out.cw 


Tin.cw)  / ln[(  T 


(33) 


hg^  = NUj.^  / dg  , 


(34) 


(35) 


and  the  subscript  cw  denotes  the  water  in  the  tankless  coil. 

3.1.3.  Boiler  Control 

The  burner  and  the  space  heating  water  circulating  pump  are  controlled  by 
on/off  control  at  upper/lower  setpoints.  A high/low  temperature  limit  switch 
for  space  heating  load  control  governs  the  circulating  pump  with  an  option 
of  manual  override.  When  the  boiler  water  temperature  was  greater  than,  or 
equal  to,  the  upper  setpoint,  the  burner  was  turned  off.  When  the  boiler 
water  temperature  was  less  than,  or  equal  to,  the  lower  setpoint,  then  the 
burner  was  turned  on.  Similarly,  when  the  space  load-side  temperature  was 
greater  than,  or  equal  to,  the  upper  limit,  the  pump  was  turned  off.  When 
the  temperature  was  less  than,  or  equal  to,  the  lower  limit,  the  pump  was 
turned  on.  In  addition,  in  order  to  achieve  better  convergence  of  computer 
simulations  with  the  HVACSIM^  program,  a very  small  amount  ( typically  0.7  % 
of  the  total  flow  rate)  of  circulating  water  was  allowed  to  flow  through  the 
pump  even  if  the  limit  controller  called  for  no  water  flow. 
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3.2.  Computer  Simulation  Procedures 
3.2.1.  Source  Prnpran 

The  boiler  component  models  described  previously  were  coded  in  the  Fortran 


77  language  as 

subroutines  bearing  names  as: 

TYPE62 

for  hot  water  boiler 

TYPE63 

for  domestic  hot  water  heating  coil 

TYPE64 

for  boiler  control 

The  subroutine  TYPE62  calls  many  routines  as  shown  in  Figure  9.  Brief 
descriptions  of  these  routines  are  given  below. 


BLINIT: 

setting  initial  conditions  at  full  load 

BLFLD: 

boiler  heat  exchanger  performance  at  full  load 

BLHX: 

boiler  heat  exchanger  performance  at  part  load 

BLFON: 

boiler  fire-box  performance  during  the  on-period  of  the 

burner 


BLOFF: 

boiler  fire-box  performance  during  the  off-period  of  the 

burner 


CPF: 

specific  heat  of  the  combustion  products 

GEF: 

gas  emissivity 

GS: 

radiation  exchange  area 

HCOF: 

convective  heat  transfer  coefficient 

PRDPP: 

mass  ratios  of  combustion  product 

PRDPR: 

viscosity  and  thermal  conductivity  of  combustion  product 

TAFF: 

adiabatic  flame  temperature 

CPCVA: 

specific  heat  of  air 

WCP: 

specific  heat  of  water 
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Figure  9 The  routines  called  by  the  subroutines,  TYPE62  and  TYPE63 
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All  routines  called  by  the  TYPE62  subroutine  are  included  in  Appendix  A, 
except  for  CPCVA  and  WCP,  which  are  part  of  the  TYPES  subroutine  in  the 
HVACSIhf'’  program. 

The  subroutine  TYPE63  needs  four  routines: 


WMU: 

viscosity  of  water. 

WK: 

thermal  conductivity  of 

water , 

WCP: 

specific  heat  of  water. 

and 

HCOVF: 

convective  heat  transfer  coefficient 

The  subroutine  TYPE64  does  not  call  any  other  routine.  The  WMU  and  WK 
routines  are  also  included  in  the  HVACSIM^  TYPES  subroutine.  The  TYPE63 
and  TYPE64  are  included  in  Appendix  A. 

3.2.2.  Input  Data  Preparation 

As  shown  in  Figure  10,  UNIT  numbers  were  assigned  to  three  component  models 
( TYPE62,  TYPE63,  and  TYPE64  ) and  index  nvimbers  were  assigned  to  state 
variables  according  to  the  HVACSIif''  program  documentation  [15,16,17].  The 
characters,  P,  M,  T,  and  C,  in  Figure  10  represent  pressure,  mass  flow  rate, 
temperature,  and  control  variables,  respectively. 

The  simulation  setup  was  accomplished  by  invoking  the  front  end  program 
HVACGEN,  which  is  included  in  the  HVACSIM^  program  package.  A simulation 
work  file  was  generated.  The  hierarchical  structure  of  the  simulation  work 
file  contains  SUPERBLOCK,  BLOCK,  and  UNIT.  In  this  boiler  simulation, 
however,  there  is  only  one  SUPERBLOCK  that  has  only  one  BLOCK  containing  three 
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Q2  = Qhwt 


Q3  = Q8ik 


^4  = Qinput 
Qs-Qlat 


^1  =^stk 
''’4  = '*’bw 
T2-T„ 
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Ts-Tnj 


Qfi-Qj 


T- 
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Figure  10  Block  diagram  of  the  boiler  model  for  use  with  HVACSIM'^ 
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UNITS.  This  simulation  setup  is  shown  in  Appendix  B.  It  was  generated  using 
the  "View  All"  command  in  HVACGEN. 

The  model  definition  file,  which  is  required  as  an  input  file  to  the  main 
program  MODSIM,  is  created  by  the  program  SLIMCON.  Whenever  information 
required  for  the  simulation  work  file  is  changed,  the  work  file  and  the 
model  definition  file  are  updated. 

A boundary  data  file,  whether  used  or  not,  must  be  provided  prior  to  a 
simulation  along  with  the  model  definition  file.  The  boundary  data  file 
may,  however,  be  empty.  For  the  simulations  in  this  study,  the  on-  and  off- 
times  of  the  water  pump  and  the  flow  rate  of  domestic  hot  water  drawn  were 
included  in  the  boundary  data  file.  Because  large  step  changes  of  a variable 
in  the  boundary  data  file  often  induce  instability  in  a simulation,  all 
large  step  changes  were  approximated  by  a number  of  small  incremental 
changes.  A program,  CRBND,  generating  such  incremental  step  changes  in  a 
boundary  data  file  is  listed  in  Appendix  C. 

The  process  of  making  a boundary  data  file  is  reasonably  simple  using  the 
program  CRBND,  when  either  space  water  heating  only  or  domestic  hot  water 
heating  only  is  being  performed,.  However,  when  space  heating  and  domestic 
water  heating  are  combined,  the  output  of  CRBND  must  be  manually  edited 
with  great  care. 

3.2.3.  Execution  of  the  Simulation 

The  main  simulation  program,  MODSIM,  was  compiled  using  an  optimized  Fortran 
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77  compiler  and  was  linked  with  necessary  library  routines  on  a mini- 
computer. The  MODSIM  program  calls  the  equation  solver  routine  that  uses 
the  Newton-Gauss  method.  With  this  method,  good  estimation  of  initial 
conditions  is  essential  to  achieve  good  convergence  of  the  solution. 

Simulations  were  performed  using  two  input  files;  the  model  definition  file 
and  boundary  data  file.  The  minimum  and  maximum  time  steps  for  simulations 
were  assigned  0.5  sec  and  200  sec,  respectively.  The  MODSIM  program 
automatically  chooses  the  time  step  and  order  of  integration  between  these 
two  limits. 

3.2.4.  Analysis  of  Simulation  Outputs 

Computer  simulation  results  were  analyzed  using  a number  of  small  programs. 
During  this  post-processing  phase,  the  simulation  output  files  were  reformatted 
to  be  usable  by  a graphics  routine,  the  heat  transfer  rates  were  integrated 
to  obtain  energy  values,  and  the  load  factors  and  fuel  utilization  efficiencies 
were  calculated.  When  interpolation  of  data  was  needed,  a routine  implementing 
a cubic  B-spline  method  was  employed.  In  addition,  a commercial  spreadsheet 
program  was  also  used  to  analyze  some  of  the  simulation  outputs. 

4.  RESULTS  AND  DISCUSSION 

Using  laboratory  test  results,  the  boiler  model  with  the  tankless  coil  was 
simulated  and  tuned.  The  tuned  model  was  then  verified  with  additional 
laboratory  measurements.  After  verification,  the  computer  model  was  used  to 
simulate  the  boiler  operation  in  both  the  summer  and  winter  seasons. 
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4. 1 Computer  Model  Verification 


Computer  simulations  for  heat-up  and  cool-down  operations  were  repeated  without 
an  external  load  until  reasonably  good  agreement  was  reached  between 
computer  simulation  results  and  experimental  measurements.  Figure  11  shows 
the  stack  gas  and  boiler  water  temperatures  of  the  computer  simulation  and  the 
laboratory  measurements  during  the  heat-up  and  cool-down  periods.  The 
integration  constants  Icfg.on*  ^cfg.off  ’-ised  equations  (26), 

(29),  and  (30)  were  determined  as  a part  of  this  process. 

Figure  12  shows  a simulation  of  part- load  operation  in  space  heating  only  mode. 
Cyclic  stack  gas  and  boiler  water  temperatures  are  compared.  In  this  case, 
space  heating  water  was  circulated  continuously  and  no  domestic  hot  water 
was  withdrawn.  The  burner  was  turned  on  or  off  depending  upon  the  boiler  water 
temperature.  It  should  be  noted  that  the  measured  boiler  water  temperature 
was  not  an  average  boiler  water  temperature  but  a local  temperature  measured 
at  the  location  of  the  temperature  sensor.  The  boiler  temperature  predicted 
by  the  computer  simulations  represented  the  average  temperature  throughout 
the  boiler.  From  the  cyclic  operation  simulations,  some  further  adjustment 
were  made  to  the  time  constants  for  the  stack  gas  and  boiler  water  during 
burner  on-  and  off -periods. 

The  simulated  pattern  of  domestic  hot  water  use  was  investigated  following 
the  test  procedure  in  the  ANSI/ASHRA.E  124P  proposed  standard  [1].  An  18- 
hour  standby  period  followed  by  a 6-hour  draw  period  were  studied  using 
laboratory  tests  and  computer  simulations.  The  stack  gas  and  boiler  water 
temperatures  are  compared  in  Figure  13.  Figure  14  depicts  the  boiler  water 
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Figure  11  The  stack  gas  and  boiler  water  temperatures  during  the  heat-up 
and  cool-down  periods 
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Temperature  (C) 


Domestic  Hot  Water  Draw  Simulation 


Figure  13  The  stack  gas  and  boiler  water  temperatures  of  a domestic  hot 
water  draw  simulation 


3A 


Boiler  Water  during  Standby 


Figure  14  The  boiler  water  tenperattire  changes  during  a part  of  the 
standby  period 


Domestic  Hot  Water  Draws 


Figure  15  The  tankless  coll  outlet  teaperattires  during  a part  of  the  draw 
period 
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Table  1 24-hour  Draw  Test 


EVENT  OIL  Qinp.exp  Qlnp.slm  ELECTRIC  Qele.exp 


(kg) 

(kJ) 

(kJ) 

(Wh) 

(kJ) 

Standby 

1.455 

65389 

67310 

78.1 

281.2 

Reheat 

0.140 

6292 

6714 

7.2 

25.9 

Draw  3 

0.500 

22471 

22800 

27.1 

97.6 

Draw  4 

0.550 

24718 

22800 

29.8 

107.3 

Draw  5 

0.540 

24268 

22800 

29.5 

106.2 

Draw  6 

0.555 

24942 

22800 

29.5 

106.2 

Draw  7 

0.545 

24493 

22800 

29.9 

107.6 

Draw  8 

0.545 

24493 

22800 

29.4 

105.8 

Total 

4.83 

217065 

210824 

260.5 

937.8 

Ratio  of 

electricity  to  oil 

— 

0.004320 
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temperature  changes  with  respect  to  time  during  a part  of  the  18 -hour  standby 
period,  while  Figure  15  shows  the  tankless  coil  output  temperatures  during  a 
part  of  the  6-hour  draw  period.  The  results  of  computer  simulations  and 
experiments  agree  reasonably  well.  Fuel  input  energy  consumption  was  used 
as  one  of  the  key  variable  in  the  comparison  of  the  laboratory  measurements 
and  the  computer  model  simulation  outputs.  Table  1 shows  fuel  consumption 
figures . 

4.2.  Simulations  of  Sitiwmer  Operating  Mode 

Using  the  verified  computer  model,  simulations  of  summer  operation  were 
performed  to  determine  input  energy  consumptions  at  various  water  loads  and 
to  evaluate  the  effect  of  boiler  room  temperature  on  the  standby  loss. 
Integration  of  power  to  obtain  energy  was  achieved  using  a utility  program, 
ENERGY.  This  program  used  a simple  trapezoidal  Integration  routine  and  is 
included  in  Appendix  D. 

Figure  16  shows  the  energy  factor  defined  as  the  ratio  of  output  energy  to 
input  energy  without  space  heading,  as  a function  of  the  domestic  hot  water 
load.  The  values  appearing  in  this  figure,  however,  were  calculated 
without  taking  into  account  the  electrical  energy  consumption.  According 
to  laboratory  measurements,  the  electrical  energy  consumed  by  the  oil  burner 
represented  approximately  0.43  % of  the  total  energy  input  for  a 24 -hour  test 
period  (see  Table  1). 

Standby  loss  variation  due  to  boiler  room  air  temperature  was  also  studied. 

A series  of  computer  simulations  were  made  with  different  room  air  temperatures. 
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Figure  16  Energy  factor  for  different  domestic  hot  water  use  level. 

Note  that  typical  average  daily  water  use  specified  by  DOE 
Test  Procedure  for  water  heaters  is  64.3  gal/day. 


38 


The  ratio,  in  percent,  of  the  energy  consumption  for  maintaining  the  boiler 
water  temperature  within  the  specified  range  during  the  standby  period  to 
the  fuel  input  energy,  Qmput , sb/Qinput*  is  shown  in  Figure  17.  The  ratio  is 
plotted  against  a dimensionless  parameter  that  is  the  difference  of  room  air 
temperature  and  boiler  water  temperature  normalized  by  the  boiler  water 
temperature,  (T^^  - T^a^/^bw*  Also  Figure  18  shows  the  ratio,  in  percent, 
of  jacket  loss  to  input  energy  against  the  same  dimensionless  parameter. 
Figures  17  and  18  show  that  the  room  air  temperature  does  affect  boiler 
standby  loss.  Laboratory  tests  were  attempted  to  confirm  these  simulation 
results,  but  no  conclusion  reached  due  to  difficulty  in  maintaining  the 
laboratory  room  air  temperature  constant  for  an  18-hour  standby  period. 

4.3.  Simulations  of  Winter  Operating  Mode 

The  effect  of  a combined  space  heating  and  domestic  hot  water  heating  load 
on  boiler  performance  was  also  studied  using  computer  simulations.  Prior  to 
the  simulations,  a typical  thermostat  cycle  rate  curve  to  meet  space  heating 
loads  was  selected  from  the  report  by  Kao,  Mastascusa,  and  Chi  [18].  The 
circulating  pump  on-  and  off-periods  were  determined  for  various  cycle  rates 
using  this  curve.  A boundary  data  file  was  generated,  which  incorporated 
the  pximp  on/off  periods  used  in  the  simulations. 

The  space  load  factor  is  given  by: 

= Qs.out  / (’?•#  Qinput)' 

where  Q,  ^ut  average  delivered  heat  flow  rate  to  the  space  during 

the  simulation  period. 
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QJ.sb/Qinput  (%)  Qi nput . sb/Qi nput  (X) 


Room  Air  Temperature  Effect 


Figure  17  The  energy  consumption  during  the  standby  period  shoving  the  room 
temperature  effect 


Figure  18  The  jacket  loss  during  the  standby  period  showing  the  room 
temperature  effect 
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Figure  19  represents  a curve  of  fuel  utilization  efficiency,  rj^  , with 
respect  to  space  load  factor,  for  space  heating  only.  This  curve  was 

created  by  smoothing  raw  simulation  outputs  of  MODSIM  by  using  the  program  SPLINE 
(see  Appendix  E) . 

With  combined  loads  of  space  heating  and  water  heating,  computer  simulations 
were  carried  out  for  12  Instead  of  24  hours  to  save  computational  time.  The 
first  6-hour  period  was  for  space  heating  only,  while  the  second  6-hour  period 
was  for  combined  space  and  water  heating.  The  energy  outputs  of  the  first 
period  were  then  multiplied  by  a factor  of  three  to  obtain  the  energy  values 
during  the  18 -hour  period  during  which  there  is  no  water  heating  load. 

Since  the  room  air  temperature  was  assumed  to  be  constant,  this  multiplication 
is  justified.  In  addition,  the  boiler  inlet  water  temperature  (returning 
from  the  heated  space)  and  the  tankless  coil  inlet  water  temperature  were 
modelled  as  remaining  constants  in  the  simulations. 

The  fuel  utilization  efficiency  for  combined  loads  is  plotted  against  the 
space  load  factor  in  Figure  20.  The  amount  of  hot  water  drawn  daily  is 
shown  as  a parameter.  At  = 0,  the  efficiency  is  the  energy  factor, 

EF.  From  this  figure,  we  can  see  that  the  water  load  causes  the  efficiency 
to  increase  when  the  space  load  is  low.  This  efficiency  increase  diminishes, 
however,  as  the  space  load  increases. 


5. 


A METHOD  FOR  D] 


NING  THE  COMBINED  SEASONAL  EFFICIENCY 


On  the  basis  of  computer  simulations  and  laboratory  experiments,  a simple 
method  for  determining  the  combined  seasonal  efficiency  of  Type  I appliances 
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Efficiency 


Interpolation 


Figure  19  Fuel  utilization  efficiency  with  respect  to  the  space  load  factor 
for  space  heating  only 


Simulation 


Spoce  Load  Factor 

Figure  20  Fuel  utillzetlon  efficiency  for  coabined  loads 
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is  presented  here.  The  primary  design  function  of  Type  I appliances  is 
space  heating  and  the  secondary  function  is  domestic  water  heating. 

Step  1:  Generate  a fuel  utilization  part-load  efficiency  curve  as  a function 
of  space  heating  loads  for  zero  domestic  water  load  having  the  form: 

17  = (a  X / (x  + b))  + c,  (37) 

and  satisfying  the  following  conditions: 

17  = 0.0  atx  = 0 

= »?225  at  X = 0.225 

17  = 17,,  atx=1.0. 

In  equation  (37),  x is  the  space  load  factor,  17225  part-load  efficiency 

at  22.5  % load  obtained  using  ANSI/ASHRAE  103  Standard,  >7,,  is  the  steady- 
state  efficiency,  and  a,  b,  and  c are  constants.  For  Type  I appliances 
these  constants  can  be  calculated  as  follows: 

a = (17..  - c)  (1  + b) 

b = -0.225  (f7,,  - 17225)  / (0-225  .7.,  - »?225  + 0.775) 
c = 0 

Figures  21  and  22  show  efficiency  curves  resulting  from  the  computer 
simulation  and  the  equation  (37)  at  60  and  120  gallons  of  daily  domestic  hot 
water  use,  respectively. 
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Ef  f I c I ency  Curve 


Figure  21  Efficiency  curves  from  computer  simulation  and  the  equation  at 
U = 60  gal 


Efficiency  Curve 


Space  Lood  Foctor 


Figure  22  Efficiency  curves  from  computer  simulation  and  the  equation  at 
U = 120  gal 
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Step  2:  Using  the  equation  (37),  find  the  efficiency,  tjI  , at  the  load 

x=0.225  + x^at«r'  where  the  average  domestic  water  load  factor,  x„ater>  given 

by: 

^^•ater  ~ Qw.out  / Qinput^*  (38) 

where  Q^.out  average  heat  flow  rate  for  domestic  water  heating. 

Step  3:  Use  the  same  form  of  equation  (37)  to  generate  a curve  of  the  form: 

fj'  = (a'x  / (X  + b'))  + c'  , (39) 

which  satisfies  three  conditions: 

= EF  at  X = 0 
17'  = »7*  at  X = 0.225 

17'  = 17,,  atx=1.0 

This  curve  represents  the  combined  efficiency  of  the  boiler  as  function  of 
space  heating  load  for  a daily  water  usage  corresponding  to  a domestic  water 
load  factor  of  x^.tar* 

Step  4:  It  is  needed  to  bin  on  equation  (37)  to  find  a correction  factor 
that  relates  the  trimmed  result  to  the  value  ri22.5  found  by  the  ASHRAE  103 
test  procedure.  Using  the  hourly  bin  data  given  in  the  report  by  Parken, 

Kelly,  and  Didlon  [19],  apply  the  bin  method  on  the  curve  represented  by 
equation  (39)  to  find  a combined  heating  seasonal  efficiency  for  the  appliance. 
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Figure  23  illustrates  Step  1 through  Step  3,  and  Table  2 shows  the  output 
of  a spreadsheet  program  for  the  boiler  with  a tankless  coil  when  the  daily 
usage  of  hot  water  is  120  gal. 

6.  CONCLUSION 

A residential,  fossil  fuel-fired,  hot  water  boiler  model  with  a tankless 
domestic  water  heating  coil  was  developed  based  upon  laboratory  tests,  and 
used  with  the  HVACSIhf'’  program.  The  computer  model  was  verified  with  some 
of  the  experimental  data  and  used  to  evaluate  the  boiler  perforaance  through 
simulations . 

Laboratory  tests  were  performed  mainly  to  obtain  parameter  values  used  in 
the  computer  model  verification.  The  computer  model  can  also  be  used  to 
simulate  other  residential  boilers  with/without  domestic  water  heating 
coils . 

A method  utilizing  a curve- fitted  equation  and  the  bin  method  is  proposed  to 
determine  the  combined  efficiency  of  a Type  I integrated  appliance  of  which 
primary  function  is  space  heating.  No  comparison  was,  however,  made  between 
this  procedure  and  draft  ASHRAE  124P  Standard  [1].  Thus,  there  is  a need 
for  additional  research  on  combined  space  and  domestic  water  heating 
appliances,  especially  Type  II  appliances  whose  primary  functions  are 
domestic  water  heating. 
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Figure  23  Illustration  of  Steps  1 through  3 of  the  NIST's  recommended 
procedure 
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Table  2 Bin  Data  Analysis 


BIN  DATA  ANALYSIS 

Region  IV 

Design  outdoor  temperature  = 5 *F 

Oversizing  factor  = 0.7 

Daily  domestic  hot  water  use  = 120  gal 

Tj  : Bin  temperature 

Nj  : Fractional  hours 

Xspace  : Space  heating  load  factor 

Xwater  : Domestic  water  heating  load  factor 

Xload  : Combined  load  factor  (=  Xspace  + Xwater  ) 

Effu  : Part- load  efficiency 

Bin  on  Xspace 

Effu  computed  by  EFFSPWT  program 


Bin  # 

Tj  (F) 

Nj 

Xspace 

Xload 

Nj*Xload  Nj*Xload/Effu 

Xwater 

Effu 

1 

62 

0.132 

0 

0.0294 

0 

0.0525 

0.0039 

0.0064 

0.0231 

0.0231 

0.4550 

0.6085 

2 

57 

0.111 

0.0784 

0.1015 

0.0087 

0.0123 

0.0231 

0.7055 

3 

52 

0.103 

0.1275 

0.1506 

0.0131 

0.0175 

0.0231 

0.7483 

4 

47 

0.093 

0.1765 

0.1996 

0.0164 

0.0213 

0.0231 

0.7723 

5 

42 

0.100 

0.2255 

0.2486 

0.0225 

0.0286 

0.0231 

0.7878 

6 

37 

0.109 

0.2745 

0.2976 

0.0299 

0.0375 

0.0231 

0.7985 

7 

32 

0.126 

0.3235 

0.3466 

0.0408 

0.0506 

0.0231 

0.8064 

8 

27 

0.087 

0.3725 

0.3956 

0.0324 

0.0399 

0.0231 

0.8125 

9 

22 

0.055 

0.4216 

0.4447 

0.0232 

0.0284 

0.0231 

0.8173 

10 

17 

0.036 

0.4706 

0.4937 

0.0169 

0.0206 

0.0231 

0.8212 

11 

12 

0.026 

0.5196 

0.5427 

0.0135 

0.0164 

0.0231 

0.8244 

12 

7 

0.013 

0.5686 

0.5917 

0.0074 

0.0089 

0.0231 

0.8271 

13 

2 

0.006 

0.6176 

0.6407 

0.0037 

0.0045 

0.0231 

0.8294 

14 

-3 

0.002 

0.6667 

0.6898 

0.0013 

0.0016 

0.0231 

0.8314 

15 

-8 

0.001 

0.7157 

0.7388 

0.0007 

0.0009 

0.0231 

0.8331 

SUM  0.234558  0.295306 

Seasonal  Efficiency  = 0.794289 
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APPENDIX  A Boiler  Component  Routines  - TYPE62,  TYPE63,  TYPE64 


SUBROUTINE  TYPE62(XIN. OUT. PAR. SAVED. lOSTAT) 


TYPE  62:  HOT  WATER  BOILER  WITH  A DOMESTIC  HOT  WATER  HEATING  COIL 


January  21.  1988  Ch«ol  Park 
Revised:  June  22.  1988 


INPUTS: 

TSTK 

TRA 

TOA 

TBLW 

PSW 

WSW 

TRW 

QCW 

SWCH 


Stack  gas  temperature  (C) 

Room  air  temperature  (C) 

Outdoar  air  temperature  (C) 

Boiler  water  temperature  (C) 

Supply  boiler  water  pressure  (kPa) 

Supply  boiler  woter  flowrate  (kg/s) 
Return  boiler  water  temperature  (C; 

Heat  flow  rate  to  the  tankless  coil  (kW) 
Switch  for  controlling  on/off 
1 for  burner  on  and  8 for  burner  off 


OUTPUTS : 

TSTK  Stock  gos  temperature  (C)  Diff.  eq. 

TBLW  Boiler  water  temperoture  (C)  Diff.  eq. 

TSW  Supply  water  temperoture  (C) 

QHWT  ■•eat  flow  to  woter  (kW) 

QSTK  Heat  loss  from  the  stack  (kW) 

OINPUT  Input  heat  flow  (kW) 

OLTNT  Latent  heat  loss  (kW) 

QJAKT  Heat  loss  from  the  boiler  jacket  (kW) 

OSS  Heat  flow  to  boiler  water  (kW) 

PARAMETERS: 


VGF  Volume  of  gas  in  the  boiler  fire-box  (m3) 

ASF  Boiler  fire-box  effective  radiation  heot  transfer  area  (m2) 
ARF  Boiler  fire-box  refractory  surfoce  orea  (m2) 

LGPF  Gas-path  length  in  the  fire-box  (m) 

EMISF  Fire-box  surface  emissivity  in  fraction 
VWB  Volume  of  woter  in  the  boiler  (mS) 

AJAKT  Boiler  jacket  surface  orea  (m2) 

UJAKT  Boiler  jocket  U-foctor  (kW/m2-C) 

CARB  Atomic  ratio  of  corbon  in  fuel 
HYOR  Atomic  ratio  of  hydrogen  in  fuel 
OXYG  Atomic  ratio  of  oxygen  in  fuel 
XNTR  Atomic  ratio  of  nitrogen  in  fuel 
SULF  Atomic  ratio  of  sulfur  in  fuel 
HFUEL  Fuel  higher  heating  value  (kJ/kg) 

CPFUEL  Fuel  specific  heot  value  (kJ/kg^) 

WFUEL  Fuel  supply  rote  (•‘fl/*) 

TFUEL  Fuel  temperature  (C) 

XAIR  Excess  air  for  combustion  (-) 

EFFYSS  Steody-stote  boiler  efficiency  in  fraction  (-) 

TBLWS  Boiler  water  temperqture  at  steady  state  (C) 

TSGS  Stock  gas  temperoture  at  full  load  (C) 

ICFGON  Integration  control  foctor  for  on-period  stock  gas  temp. 
ICFGOF  Integrotion  control  foctor  for  off-period  stock  gas  temp. 
ICFW  Integration  control  factor  for  water  temperoture 
DTFG  Short  stort  period  for  burner  on/off  cycle  (s) 

ICFGNS  Integration  control  for  short  start  on-period 
ICFGFS  Integration  control  for  short  stort  off-period 
OTBW  Minimum  time  interval  for  updating  QSS  (s) 


C 


PARAMETER  (NSAVED-5 .NDE-2 .NIN-9 .NOUT-9 .NPAR-28) 

LOGICAL  COMINT 

REAL  LGPF . ICFGON . ICFGOF . ICFW. ICFGNS . ICFGFS 

DIMENSION  XIN(NIN) .OUT(NOUT) .PAR(NPAR) . lOSTAT(NOUT) . SAVED (NSAVED) 
COMyON  /CHRONO/  TIME.TSTEP.TTIME.TMIN. ITIME 
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COtAdON 

COMMON 

COMMON 

COMMON 

DATA 

DATA 

C»  DATA 

C»  NAMELIST 

C*  NAMELIST 

C»  4 
C*  NAMELIST 

C*  NAMELIST 

C*  NAMELIST 

C*  NAMELIST 


/PRODCT/  PT , PC02 . PH20 , PN2 . P02 . PS02 . RATE . RWTF , RPTF 
/CONFIG/  VGF.ASF.ARF, LGPF.EMISF.WVB.AJAKT.UJAKT 
/FUEL/  HFUEL.CPFUEL.WFUEL.TFUEL 
/FULLD/  CNTU.TSGS.QINP.TAUG.TAUW 

COMINT/. TRUE./ 

WTFAC  /0.8/.TBLREF/97.5/ 

WTFAC  /0 . 8/ . TRAREF/23 . 0/ . TBLREF/97 . 5/ 

/NAM1/  OINP 

/NAM2/  TIME. SWCH . OSS .TSTK.TBLW.OF. OHXC , 

DTSTK . DTB  LW . TEXGF . WGOFF . WTT  L 
/NAM3/  TIME. IONOFF.TIME0,DELT,ICFGON. ICFGOF 
/NAM4/  TIME. TIME1 .DELTB. OSS. OSS1 
/NAM5/  TIME.TBLW.TSTK.OCW 

/NAM6/  T I ME . 01 NPUT . OLTNT . OSTK . Q J ART . QHWT . OSS 


I nput  s ; 

TSTK  =XIN(1) 

TRA  =XIN(2) 

TOA  =XIN(3) 

TBLW  =XIN(4) 

PSW  =XIN(5) 

WSW  =XIN(6) 

TRW  =XIN(7) 

OCW  =XIN(8) 

SWCH  =XIN(9) 

Poromete  rs : 

VGF  -PAR(1) 

ASF  -par (2) 

ARF  =PAR(3) 

LGPF  =PARf4) 

EMISF  =PAR(5) 

VWB  -PAR(6) 

AJAKT  -PAR(7) 

UJAKT  -par (8) 

CARB  =PAR(9) 

HYDR  =PAR(10) 

OXYG  =PAR(1l) 

XNTR  =PAR(12) 

SULF  =PAR(13) 

HFUEL  -PAR(14) 

CPFUEL  -PAR(15) 

WFUEL  -PAR(16) 

TFUEL  =°AR(17) 

XAIR  =PAR(18) 

EFFYSS  -=PAR(19) 

TBLWS  -PAR(20) 

TSGS  -FAR (21) 

ICFGON  -PARf22) 

ICFGOF  -PAR (23) 

ICFW  -PAR(24) 

DTFG  -PAR(25) 

ICFGNS  -PAR (26) 

ICFGFS  -PAR (27) 

DTBW  -PAR (28) 

Assume  that  the  boiler  supply  water  temperature  is  the  same 
as  the  boiler  water  temperature,  i.e.  TSW-TBLW. 

TSW-TBLW 


Initial  conditions  ot  steady  state.  Input  data  at 
beginning  of  simulation  must  be  entered  using  o steady— stote 
va I ues . 

1F(ITIME.EQ. 1 .AND.  COMINT)  THEN 
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c* 


PT=1  . 

TSTK=TSGS 

TBLW=TBLWS 

TSWS=TBLWS 

CALL  BL I NI T (CARB , HYDR . OXYG . XNTR , SULF , XAIR , 
k TRA.TBLWS.TSWS.EFFYSS.WGON) 

COM  I NT=. FALSE. 

SAVED(1 )=TIME 
SAVED(2)=WG0N 
SAVED(3)=0.e 
SAVED(4)=1 
PRINT  NAM1 

ELSE 

WG0N=SAVED(2) 

END  IF 


Burner  off 

IF(TIME  .GT.  SAVED(I))  THEN 
IF(SWCH, LE.0.5)  then 

CALL  BLFOFF(TSTK,TRA.TBLW,VyGON.WGOFF.QFC.TEXGF) 
CALL  BLMX(TSW.TEXGF,TSTK.WGOFF.OHXC) 

CALL  CPCVA(TRA,CPA.DUI^,DUMMY,DUIvWY) 

OF=OFC 

IONOFF*0 

CAPG=WG0FF*CPA 


Burner  on 
ELSE 

CALL  BLFON( TRA . TB LW , WGON , OFCR . TAFC , TEXGF) 
CALL  BLHX(TSW. TEXGF. TSTK, WGON, OHXC) 
QF=QFCR 
IONOFF«1 

CAPG=WGON*CPF(TRA.TSTK) 

END  IF 

END  IF 


Heat  fluxes 

CAPW=WSW«WCP(TBLW) 

OlNPUT=WFUEL*HFUEL*IONOFF 
OLTNT=2442 . •WFUEL*RWTF* lONOFF 
0 J AK  T=A J AKT • U J AKT • ( TB  LW-TRA ) 
QHWT=CAPW* ( TSW-TRW)+OCW 
CSTK=QINPUT-OLTNT-OF-OHXC 
oss=qf+ohxc-qjakt-ocw 
OSS1=SAVED(3) 


Setting  up  a short  time  intervol  of  rapid  change  of  gas 
temperature  end  using  different  integration  constants 

IF(TIME  .GT  SAVED(1))  THEN 
I BURN=SAVED ( 4 ) +0 . 00 1 
1F(I0N0FF  .NE.  IBURN)  THEN 
TIME0  - TIME 
ELSE 

DELT  - TIME-TIME© 

END  IF 

SAVED(4)  - lONOFF 
END  IF 

IF(DELT  .LT.  DTFG)  THEN 
ICFGON  - ICFGNS 
ICFGOF  « ICFGFS 
END  IF 


Introducing  time-deloy  to  the  boiler  water  temperature. 
OSS  is  updated  every  t ime— i nte rva I of  DELTB,  if  DELTB  is 
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greater  than  ar  equal  to  the  specified  minimum  value  of  DTBW. 

At  the  some  time,  the  reference  time.TIMEI, 

is  olso  changed.  It  is  assumed  that  OSS  holds  a steady  value 
in  the  t i me— i n te rvo I . 

NOTE  thot  DELTB  depends  upon  the  specified  maximum  time  step, 

TMAX . Due  to  this  foct,  the  use  of  different  TMAX ’ s could  result 
in  different  simulotion  results. 

1F(1TIME  .EQ.  1)  THEN 
TIME1=TIME 
SAVED(5'=T1ME1 
END  IF 

IF(TIME  .GT.  SAVED(I))  THEN 
TIME1=SAVED(5) 

DELTB  = TIME-TIME1 
IF(DELTB  .GE.  DTBW)  THEN 
SAVED(3)=OSS 
TIME1=TIME 
SAVED(5)=TIME1 
END  IF 
END  IF 


Derivatives  of  stock  gos  ond  boiler  water  temperotures 

IFdONOFF  .EQ.  1)  THEN 

DTSTK=(TSGS-TSTK)/(ICFGON.TAUG) 

ELSEIF(IONOFF  .EO.  0)  THEN 

T I N F=WT  FAC • TB  LR  EF+ ( 1 . -WT  FAC ) • TRA 
DTSTK*(TINF-TSTK)/( ICFGOF.TAUG) 

END  IF 

DTBLW=(CAPW*(TRW-TSW)+0SS1 )/( ICF¥(*TAUW) 

Outputs 

OUT(1)  =.OTSTK 
OUT (2)  =DTBLW 
OUT (3)  -TSW 
OUT (4)  =QHWT 
OUT (5)  =QSTK 
OUT(6)  =QINPUT 
OUT(7)  «<}LTNT 
OUT (8)  =QJAKT 
OUT(9)  =OSS 

IF(TIME.  GT.  SAVED(I))  THEN 
PRINT  NAM2 
PRINT  NAM3 
PRINT  NAM4 
PRINT  NAM5 
PRINT  NAM6 
END  IF 

SAVED(1)=*TIME 

lOSTAT(1)-0 

IOSTAT(2)-0 

IOSTAT(3)»0 

I0STAT(4)-1 

IOSTAT(5)-1 

IOSTAT(6)-1 

IOSTAT(7)-1 

I0STAT(8)=1 

IOSTAT(9)-1 

RETURN 

END 


SUBROUT  I NE  B L I N I T ( CARB , HYDR . OXYG . XNTR . SU  LF , XA I R , 


54 


oooooo  oooooooo  ooo  oo  o ooo  ooo  ooooo  oooo  ooooooooo 


4 TRA.TBLW.TSW.EFFYSS.WGON) 


CNTU  Modified  boiler  HX  heat  tronsfer  number  at  full  load 

TAUG  Modified  boiler  stock  gos  time  constant 

TAUW  Modified  boiler  water  time  constant 


REAL  LGPF 

COKWON  /PRODCT/  PT . PC02 . PH20. PN2 , P02 . PS02 . RATF . RWTF . RPTF 

COMMON  /CONFIG/  VGF  . ASF  . ARF  . LGPi^  . EMISF  . VWB  , AJAKT  . UJAKT 

COMMON  /^UEL/  HFUEL.CPFUEL.WFUEL.TFUEL 

COMMON  /FULLD/  CNTU , TSGS . QINP . TAUG . TAUW 

namelist  /NAMINI/  WGON.OFCR.TAFC.TEXGF.QJAKT.QFNS. TAUG. TAUW 

Combustion  gos  properties 

CALL  PRDPP(CARB . HVDR . OXYG , XNTR . SULF , XAIR) 

Estimote  the  gos  polh  lenth  in  the  boiler  fire-box  by  matching 
the  calculated  ond  the  meosured  gas  temperotures  at  the  exit  of 
the  f i re— box . 

CALL  BLFON(TRA.TBLW,WGON.OFCR.TAFC.TEXGF) 

Steody-stote  condition 

0JAK’’=AJAKT*UJAKT*(TBLW-TRA) 

TINHX=TEXGF 

ofns=qfcr-qjakt 

CALL  B LFLD (WGON . QFNS . T INHX . EFFYSS . TSW) 

Time  constants  of  gas  and  water 

IF(TAUG.LT.1 .E-10)  THEN 
T1G=(TFUEL+RATF*TRA)/RPTF 

T1G=(CPFUEL*TFUEL+RATF»TRA)/RPTF 

T2G=TSGS 

CPG=CPF(T1G.T2G) 

CAPG=CPG«WG0N 
TAUG=1 .0/CAPG 
END  IF 

IF(TAUW.LE.1 .E-10)  THEN 
TAUW=4200. ‘MWE 
END  IF 

PRINT  NA:>(INI 

RETURN 

END 


SUBROUT  I NE  B LF LD ( WGON , OFNS , T I NHX . EFFYSS .TSW) 


Boiler  heot  exchonger  performonce  ot  full  lood.  Colculotion 
of  the  number  of  heot  transfer  unit. 


REAL  K.MU.NTU 

COMMON  /FULLD/  CNTU. TSGS. QINP. TAUG. TAUW 

NAMELIST  /NAMFLD/  OHXSS .NTU .CNTU 

The  specific  heat  of  gos  poss i ng  through  the  heot 
exchanger  by  using  the  meosured  stock  gas  temperoture  ot 
o steady  state. 

CP»CPF(TINHX.TSGS) 
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C The  number  of  tronsfer  units  at  full  lood 

C 

OHXSS=EFrYSS*OINP-C)FNS 

ENTU=1 .-QHXSS/(CP*WGON.(TINHX-TSW)) 

IF(ENTU.GT. 1 .E-20)  THEN 
NTU=LOG(1 ./ENTU) 

ELSE 
NTU=50 . 

CP=CPF(TINHX.TSW) 

TAVE=0 . 5* ( T 1 NHX+TSW) 

PRINT  TOO  LARGE  EFFICIENCY  VALUE  ’ 

END  IF 

CALL  PRDPR(TAVE.MU.K) 
CNTU=NTU*WGON«*0.2*MU**0.4*(CP/K)»*0.6 
C«  PRINT  NAMFLD 

C 

RETURN 

END 

C 

C*  


SUBROUT  I NE  B LHX  ( T Svy , T I NHX  . T EXHX . WG  . GHXC ) 

C 

c 

c 

C Boiler  heat-exchcnger  port  lood  performance 

C 



c 

REAL  K.’.fU.NTU 

COMMON  /FULLD/  CNTU . TSGS .OINP , TAUG . TAUW 

C 

C The  number  of  tronsfer  unit  for  boiler  heot  exchonger 
C at  port  load  condition. 

C 

TAVE=0 . 5* (T I NHX+TEXHX) 

1F(ABS(TINHX-TEXHX)  .GT.  0.001)  THEN 
CP=CPF(TINHX.TEXHX) 

ELSE 

CALL  CPCVA(TAVE. CPA. DUMMY, DUWWY.DUWWY) 

END  IF 
C 

CALL  PRDPR(TAVE.MU.K) 

C 

C Convective  heot  transfer  rote 
C 

IF(WG.GT.  1.0E-6)  THEN 

NTU=CNTU/(WG.»0.2»MU**0.4«(CP/K)**0.6) 

QHXC=CP*VVG«  ( T INHX-TSW)  • ( 1 . -EXP(-NTU) ) 

ELSE 

QHXC=CP»WG* ( T I NHX-TSW) 

END  IF 
C 

RETURN 

END 

C 



C 

SUBROUTINE  BLFON(TRA , TBLW . WGON .QFCR . TAFC . TEXGF) 

C 

C 

c 

C Simulation  of  boiler  fire-box  performance  during  on-period 
C 


c 

PT 

Gas  pressure  (otm) 

c 

PC02 

Number  of  moles 

of 

C02 

c 

PH20 

Number  of  moles 

of 

H20 

c 

PS02 

Number  of  moles 

of 

S02 

c 

P02 

Number  of  moles 

of 

02 

c 

PN2 

Number  of  moles 

of 

N2 

c 

RATF 

Moss  rat i o of  a i 

i r 

to  fuel 

c 

RWTF 

Moss  ratio  of  water  to  fue 
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RPTF  Mass  ratio  of  combustion  product  to  fuel 
TS  Fire-box  surface  temperature  (C) 

OFCR  Fire-box  heat  transfer  rate  (kW) 

TAFC  Adiabatic  f I ome  tem.perature  (C) 

TEXGF  Fire-box  exit  gas  temperature  (C) 


REAL  K.MU.LGPF 

COMMON  /PRODCT/  PT . PC02 . PH20 . PN2 , P02 . PS02 . RATF , RWTF . RPTF 
COMMON  /CONFIG/  VGF , ASF . ARF , LGPF . EMISF , VWB , AJAKT . UJAKT 
COtvtytON  /FUEL/  HFUEL.CPFUEL.WFUEL.TFUEL 

CCMJON  /FULLD/  CNTU . TSGS , QINP . TAUG , TAUW 

DATA  SIGMA/5.670E-11/.  CKELVN/273 . 1 5/ , PI/3.14159/ 


Fuel  input 

hhv=hfuel+cpfuel*(tfuel-tra) 

QINP=WFUEL*HHV 

Set  the  fire-box  surfoce  tmperature  to  be  the  same  os  the  boiler 
water  temperoture. 

TS=TBLW 


Convert  temperature  unit  into  absolute  unit 

TKS=TS+CKELVN 

TKA=TRA+CKELVN 

The  gas  and  water  vopor  flow  rotes  during  on-cycle 

WGON=WFUEL*RPTF 

WH20--WFUEL*RWTF 

HF=WFUEL*HHV-2442 . •WH20 
ATOTAL=ASF+ARF 
CS=ASF/ATOTAL 
AL=3.5»VGF/AT0TAL 

Adiobotic  f I ome  temperature 

TAFC=TAFF(HHV.TRA) 

TAF=TAFC+CKELVN 

The  heat  transfer  rote  due  to  radiotion  and 
convection  using  the  well-stirred  furnoce  theory 

TG2=TAF-250 

ITR-0 

ITR-ITR+1 

TG2C=TG2-CKELVN 

TGAVE=0.5*{TAF+TG2) 

CP0E=CPF(TRA.TG2C) 

cpae=:pf(tg2C.tafc) 

CP0A=CPF(TRA.TAFC) 

TT»TGAVE-CKELVN 
CALL  PRDPR(TT,MU.K) 

PR=MU*CPAEA 

ac=vgf/lgpf 

DH=SQRT(4.0.AC/PI) 

R E=WGON • DH/ ( MU • AC ) 

HCOV=HCOVF ( RE . PR . K , DH) 

GE=GEF(PC02,PH20.AL,TT) 

AGS=GS ( ATOT AL . CS , GE . EM I SF ) 

TRATIO=TKS/TGAVE 

AGS=AGS*(1 .-TRATIO**3)/(1 .-TRATI0*»4) 

AGS=AGS+2 . • ASF . HCOV/ ( S I GMA • ( TGAVE+TKS ) • • 3 ) 
0FCR-AGS*SIGMA.(TGAVE**4-TKS**4) 

A new  estifflote  of  gos  temperature  by  using  Newton's  method 
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CAPWOE=CPOE*WGON 
CAPWOA=CPOA*WGON 
CCOAT=CAPWOA. (TAF-TKA) . (HF-QFCR) 
CCOE=CAPWOE»HF 

tgcal=ccoat/ccoe+tka 
DTG2=TGCAL-TG2 
IF(ABS(DTG2) .GT.  10. ) THEN 
IF( ITR.GT. 1 ) THEN 

TG=TG1-DTG1»(TG2-TG1 )/(DTG2-DTG1  ) 
IF(ABS(DTG1 ) .GT.ABS(DTG2))  THEN 
DTG1=DTG2 
TG1=TG2 
END  IF 
TG2*TG 
ELSE 

DTG1=DTG2 

TG1=TG2 

TG2=TG2+100. 

Ih  fOTGI . LT.0, ) THEN 
■ 2=TG2-200. 

END  IF 
END  IF 
GOTO  10 
END  IF 


Gos  temperature  ot  boiler  fire-box  during  on-period 

TEXGF=TG2-CKELVN 

RETURN 

END 


SUBROUT  I NE  BLFOF  F ( TSTK  , TRA , TB  LW . WGON , VKGOFF , QFC . TEXGF) 


Boiler  fire-box  performance  during  off-period 


REAL  K.MU.NTU. LGPF 

COMMON  /CONFIG/  VGF , ASF . ARF . LGPF . EMI SF . VWB . AJAKT . UJAKT 
COMMON  /FULLD/  CNTU . TSGS . QINP . TAUG . TAUW 

DATA  CKELVN/273.15/.  PI/3.14159/ 

Moss  flow  rate  of  air  during  off-period 

IF(TSTK-TRA  .GT.  0.01)  THEN 

IF(ABS(TSTK-TRA)  .GT.  0.01)  THEN 
WGOFF=WGON* ( ( TSTK-TRA)/( TSGS-TRA) ) • .0 . 56 
k •((TSGS+CKELVN)/(TSTK+CKELVN))..1 .19 

DS»0.4 

WGOFF»OS»WGOFF 

ELSE 

WGOFF-0 . 0 
END  IF 


The  properties  of  gos  in  the  boiler  fire-box: 

The  dynomic  viscosity,  thermal  conductivity,  specific  heot 
copocity,  Prant I number,  ond  Reynolds  number  bosed  on  the 
hydrau lie  d i ometer . 

TAVE=0.5*(.TRA+TEXGF) 

CALL  PRDPR(TAVE,MU,K) 

CALL  CPCVA ( T EXGF , CPA , DUNWY , DUVI-ilY , DUK-WY) 

PR=MU*CPA/K 
AC=VGF/LGPF 
DH=SORT( 4.0. AC/PI) 

RE=WGOFF*DH/(MU.AC) 

The  convective  heot  tronsfer  coefficient. 
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c 


HCOV=HCOVF ( RE , PR . K . DH ) 


The  off-cycle  heat  transfer  -rote  using  the 
effectiveness  method  of  o heat  exchonger. 

The  gas  temperoture  at  the  exit  of  the  boiler  fire-box. 

CAPA=CPA*WGOFF 
IF(WGOFF  .GT.  0.001)  THEN 
NT  U=AS  F • HCOV/C AP A 

OFC=CAPA* ( TRA-TBLW) • ( 1 . -EXP (-NTU ) ) 

TEXGF=TRA-QFC*(1 ./CAPA) 

ELSE 
QFC=0 . 0 

TEXGF=(TRA+TBLW)/2. 

END  IF 

RETURN 

END 


SUBROUTINE  TYPE63(XIN.0UT,PAR.SAVED. lOSTAT) 


TYPE  63;  HOT  WATER  COIL  WITH  CONSTANT  WALL  TEMPERATURE 

February  16,  1988  Cheol  F rk 
Revised  : April  15,  1988 

INPUTS: 

TBLW  Boiler  woter  temperoture  (C) 

WCW  Woter  flow  rote  through  the  coil  (kg/e) 

TICW  Inlet  coil  woter  temperoture  (C) 

OUTPUTS: 

TOCW  Outlet  coil  woter  temperoture  (C) 

QCW  Heot  flow  rote  from  boiler  woter  to  coil  woter  (kW) 
PARAMETERS: 

DCOI L Diometer  of  coil 
LCC  L Length  of  coi I 


parameter  ( NSAVED-2 , NDE=0 . N I N=3 , NOUT=2 , NPAR*2 ) 

REAL  K,LCOIL,LMTD.MU 

DIMENSION  XIN(NIN) ,OUT(NOUT) ,PAR(NPAR) , lOSTAT(NOUT) ,SAVED(NSAVED) 
COMMON  /CHRONO/  T IME , TSTEP , TT IME , TMIN , I T IME 
DATA  PI/3.14159/ 

NAMELIST  /NAM1/  HCOV , TOCW, LMTD . OCW 
NAMELIST  /NAM2/  TIME, TBLW, WCW, TICW, TOCW, OCW 

Inputs: 

TBLW  -XIN(1) 

WCW  -XIN(2) 

TICW  =XIN(3) 

Porometers: 

DCOIL  -PAR(1) 

LCOIL  =PAR(2) 

Temporolly  set  the  wall  temperature  of  the  woter  coil  equal  to 
the  boiler  woter  temperoture. 

TS-TBLW 

Woter  properities;  dynomic  viscosity,  thermol  conductivity,  ond 
specific  heot  of  water  from  the  subroutine  WATPR  with  overage 
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C temperoture  except  at  the  beginning. 

C 

ir(ITlME.EO.I)  THEN 
T AVE=0 . 5» ( TB LW+T I CW ) 

SAVED( 1 )=T1ME 
SAVED(2)=TICW 
ELSE 

TOCW=SAVED(2) 

TAVE=0.5*(TOCW+TICW) 

END  IF 

MU=WMU(TAVE) 

K=WK(TAVE) 

CP=WCP(TAVE) 

C 

C The  coil  surface  area,  and  wetted  c ross-sec t i na I orea. 

C 

as=p:*dcoil«lcoil 

AC=P I . DCO I L • DCO I L/4 . 0 
C 

C Reynoics  number  and  Pronditl  number  of  the  coil  woter 

C 

RE=WCW.DCOlL/(AC»MU) 

PR=MU»CP/K 

HC0V=HC0VF ( re . PR . K , DCO I L ) 

C 

C The  outlet  temperature  and  the  heat  flow  rote  when 
C the  coil  surfoce  temperoture  is  constant. 

C 

IF(WCV».GT.  0.0001)  THEN 
A=HCOV.AS/(WCW*CP) 

TOCW=TS- ( TS-T I CW ) . EXP ( -A ) 

LMTD= ( TOCW-T I CW , /log ( ( TS-T I CW) / ( TS-TOCW ) ) 

OCW=HCOV*AS*LMTD 
ELSE 
TOCW=TS 
OCW-0.0 
END  IF 
C 

C Outputs: 

C 

OUT(1)  =TOCW 
OUT ( 2 ) =OCW 
SAVED(2)=TOCW 
C 

IF(TIME  .GT.  SAVED(I))  THEN 
C»  PRINT  NAM1 

C«  PRINT  NAM2 

ENDI'" 

SAVED(1)=TIME 

C 

10STAT(1 )-0 
lOSTAT(2)-0 
C 

RETURN 

END 

C 

C... 

c 

SUBROUT  I NE  T YPE64 ( X I N . OUT , PAR , SAVED . I OSTAT ) 

C 

C 

C 

C TYPE  64;  BOILER  BURNER  AND  CIRCULATING  PUMP  CONTROLS 
C 

C Febroory  17,  1988  Cheol  Pork 
C Revised  : June  6,  1988 

C 

C INPUTS; 

C TBURN  Burner  control  temperature  (C) 

C TBLW  Boiler  woter  temperoture  (C) 

C "^LOAD  Lood  temperoture  (C) 

C CTHERM  Thermostat  control  indicator  for  the  burner 

C 1 for  thermostot  control,  0 for  monuol  control  (-) 
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WCW  Domestic  hot  water  flow  rote  (kg/s) 

OUTPUTS; 

SWCH  Burner  control  on/off  signal  (-) 

WSW  Boiler  circuloting  woter  flow  rote  (kg/s) 

PARAMETERS: 

TBURON  Boiler  burner-on  temperoture  (C) 

TBUROF  Boiler  burner-off  temperature  (C) 

TPMPON  Boiler  water  circulating  pump-on  temperature  (C) 
TPMPOF  Boiler  woter  circulating  pump— off  temperature  (C) 
WPUMP  Circuloting  wcter  flow  rate  (kg/s) 

RATEOF  Minimum  woter  flow  rote  ratio 

during  off-period  of  pump  (-) 


PARAMETER  (NSAVED=1 . NDE=0 . NIN=5 , NOUT=2 . NPAR=6) 

DIMENSION  XIN(NIN) ,OUT(NOUT) .PAR(NPAR) . IOSTAT(NOUT) . SAVED (NSAVED) 
COMMON  /CHRONO/  T IME . TSTEP . TT IME . TMIN , I T IME 

NAMELIST  /NAM1/  T IME . TBUROF . TBURON . TBLW , SWCH .WCW 

1 nput  s ; 

TBURN  =XIN(1) 

TBLW  =XIN(2) 

TLOAD  =XIN(3) 

CTHERM  =XIN(4) 

WCW  =XIN(5) 

Porometers: 

TBURON  =PAR(1) 

TBUROF  =PAR(2) 

TPMPON  =PAR(3) 

TPMPOF  =PAR(4) 

WPUMP  =PAR(5) 

RATEOF  =PAR(6) 

1F(1TIME  .EQ.  1)  THEN 
SAVED(1)=TIME 
END  IF 

Boiler  burner  control 

The  burner  is  controlled  by  the  thermostot  sensing  the  boiler 
woter  temperature  if  the  indicator,  CTHERM,  is  greater  than  or 
equol  to  0.5.  Otherwise  TBURN  controls  the  burner. 

The  VO  lues  of  CTHERM  and  TBURN  must  be  present  in  the  boundory 
data  file. 

1F(CTHERM  .GE.  0.5)  THEN 
TBURN-TBLW 
END  IF 

Turn  on  ond  off  the  burner  ond  the  pump  only  when  the  time  step 
vories  to  prevent  folse  action  due  to  numerical  instobility 
during  the  lorge  chonge  of  o vorioble. 

IF(TIME  .GT.  SAVED(I))  THEN 
I F( TBURN  .GE.  TBUROF)  THEN 
SWCH*0 . 0 

ELSEIF(TBURN  .LE.  TBURON)  THEN 
SWCH-1.0 
END  IF 

C Boiler  water  circulating  pump  control 

IF(TLOAD  .LE.  TPMPON)  THEN 
WSW*WPLWP 

ELSEIF(TLOAD  .GE.  TPMPOF)  THEN 
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WSW=RATEOF»WPUMP 
END  IF 
END  IF 


Outputs ; 

OUT(1)  =SWCH 
OUT (2)  =WSW 
C 

IOSTAT(1 )=0 
IOSTAT(2)=0 
C 

C*  IF(T1ME  .GT.  SAVED(I))  THEN 

C*  PRINT  NAM1 

C»  ENDIF 

SAVED(1)=TIME 

C 

RETURN 

END 
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FUNCTION  CPF (TCI .TC2) 


Combustion  product  specific  heat 
INPUTS: 

TCI  Air  temperature  (C) 

TC2  Combustion  product  temperature  (C) 

PC02  Number  of  moles  of  C02 

PH20  Number  of  moles  of  H20 

P02  Number  of  moles  of  02 

PS02  Number  of  moles  of  S02 

PN2  Number  of  moles  of  N2 

OUTPUT : 

CPF  Specific  heot  of  combustion  product  (kJ/(kg— C)) 


COkiMON/PRODCT/PT  . PC02 . PH20 , PN2  . P02  . PS02 . RATF  , RWTF  . RPTF 

DIMENSION  C(6.2) 

DIMENSION  A(30) .0(30) 

data  A/4  2497678 , -6 . 91 2652E-3 .3.1 6021 3AE-5 . -2 . 971 5432E-8 . 
*9.510358E-12,2. 1701 . 1 . 03781 1 5E-2 . -1 . 0733938E-5 . 6 . 34591 75E-9 ,-1 .628 
t070lE-1 2. 4. 1565016,-1 . 7244334E-3. 5 . €98231 6E-6. -4. 59300448-9. 1 .4233 
ft654E-1 2 .3.6916148.-1. 3332552E-3 . 2 . 65031 E-6 ,-9 . 7688348-10,-9 . 977223 
*4E-1 4,3.7189946.-2.51 67288E-3 . 8 . 5837353E-6 , -8 . 299871 6E-9 , 2 . 7082 1 8E 
4:-1 2 ,3 . 2257132 . 5 . 6551 207E-3 .-2 . 49702088-7,-4 . 22067668-9 .2.1 3927338- 
412/ 

DATA  B/1 .1795744.1 .09505948-2,-4.0622138-6.7.13702818-10, 

4-4 . 74903538-1 4 , 4 .4129266 .3.1 9228968-3,-1 . 2978238-6 ,2.41 474468-10 
4.1. 67429868-1 4.2. 6707532 .3.031711 58-3 , -8 . 5351578-7 .1.1 7908538-1 0 . 
4-6 . 1 9735688-15 . 2 . 8545761 . 1 . 59763168-3,-6 . 25662548-7 , 1 . 13158498- 
410,-7. 6897078-1 5.3. 59761 29,7.81 456038-4,-2 . 2386678-7 . 4 . 2490 1 598 
4-1 1 . -3 . 34602048- 15.5.1982451.2. 05950958-3 . -8 . 62544508-7 . 

41 .66365238-10,-1 .18478378-14/ 

Set  T1  ond  T2  in  degree  K 

T1-273. 15+TC1 
T2-273. 15+TC2 
TM-1000. 

00  10  1-1 .6 
DO  10  J-1 .2 
C(1 .J)-0. 

10  CONTINUE 
C 

IF(T1 . LT. 1000. ) THEN 
IF(T2.LT. 1000. ) THEN 
DO  30  1-1 .6 
►4-5*(I-1) 

DO  20  J-1 .5 

C(1 .1)^(1. 1)+A(JJ)*(T2*»J-T1.*J)/J 
20  CONTINUE 

C( 1 . 1 )-C( 1 . 1 ) *8 . 32066/(T2-T1 ) 

30  CONTINUE 

CPF-(PC02*C(2 . 1 )+PH20*C(3. 1 )+PN2*C(4. 1 ) 

4 +P02*C(5. 1 )+PS02*C(6.1))/(PC02*44. 

4 +PH20* 1 8 . +PN2 • 28 . +P02 • 32 . +PS02 • 64 . ) 

ELSE 

DO  50  1-1 .6 
M-5«(I-1) 

DO  40  J-1 .5 
JJ-M+J 

C(I .1)-C(I .1)+A(JJ).(TM.*J-T1..J)/J 
C(I .2)-C(I .2)+B(JJ)*(T2**J-TM.*J)/J 
40  CONTINUE 

C(I.1)-(C(I.1)+C(I,2))*8.32066/(T2-T1) 
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CONTINUE 

CPF=(PC02»C(2. 1 )+PH20»C(3, 1 )+PN2*C(4. 1 ) 
4 +P02*C(5. 1 )+PS02*C(6. 1 ))/(PC02»44. 

4 +PH20»18.+PN2»28.+P02»32.+PS02*64. ) 

ENDIF 
ELSE 

IF(T2.GT. 1000. ) THEN 
DO  70  1=1 .6 
M=5*(I-1 ) 

DO  60  J=1 ,5 
JJ=M+J 

C(I.2)=C(I .2)+B(JJ)*(T2**J-T1*.J)/J 
60  CONTINUE 

C(I,2)=C(I ,2)*8.32066/(T2-T1) 

70  CONTINUE 

CPF=(PC02«C(2.2)+PH20*C(3.2)+PN2*C(4,2) 
4 +P02«C(5.2)+PS02.C(6.2))/(PC02»44. 

4 +PH20*18.+PN2*28.+P02*32.+PS02*64. ) 

ELSE 

DO  90  1=1 ,6 
W=5«(I-1) 

DO  80  J=1 .5 
J J=VM-J 

C(I . 1 )=C( I ,1 )+A(JJ)*(TM..J-T2**J)/J 
C(I .2)=C(1 .2)+B(JJ)*(T1..J-TM».J)/J 
80  CONTINUE 

C(I.1)=(C{I.1)+C(I .2))»8.32066/(T1-T2) 
90  CONTINUE 

CPF=(PC02«C(2. 1)+PH20*C(3. 1)+PN2*C(4, 1) 
4 +P02 • C ( 5 . 1 ) +PS02  *0 ( 6 . 1 ) ) / ( PC02 • 44 . 

4 +PH20*18.+PN2*28.+P02*32.+PS02»64. ) 

ENDIF 
ENDIF 
RETURN 
END 


FUNCTION  GEF(PPC02.PPH20.XL.TG) 


Gas  emissivity 
INPUTS: 

PPC02  Portiol  pressure  of  C02  (otm) 

PPH20  Portiol  pressure  of  H20  (otm) 

XL  Length  of  combustion  (ft) 

TG  Rodioting  gos  temperoture  (C) 

OUTPUT : 

GEF  Gos  emi ss i V i ty 


Select  gos  emissivity  constants  C1  4 C2 

XX-3 . 2808*XL* (PPC02+PPH20) 

IF(XX.GT.0.2)  THEN 
C1-287. 

C2=0.4 

ELSE 

C1=406.9 

C2=0.62 

ENDIF 

Compute  gos  emissivity 

TRG-1 .8*(TG+273. 15) 

GEX«XX**C2«C1/TRG 

GEF-GEX 

RETURN 

END 
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FUNCTION  GS(AT. CS.GE.se) 


Gos/Sink  exchange  area 
INPUTS; 

AT  Total  areo  (tn2) 

CS  Ratio  of  active  furnoce  surface  areo  to  total 
GE  Gas  etn  i S3  i V i ty 

SE  Active  surface  emissivity 

OUTPUT : 

GS  Rodiotion  exchange  oreo  (m2) 


area 


GSX=1 ./GE+1 ./(CS*SE)-1 . 

GSX=AT/GSX 

GS=GSX 

RETURN 

END 


FUNCTION  HCOVF(RE.PR.K.D) 


Convective  heot  tronsfer  coefficient 
INPUTS: 

RE  Reynolds  number 

PR  Prondt I number 

K Thermol  conductivity  (kW/(m— C)) 

D Hydoul ic  diometer  (m) 

OUTPUT : 

HCOVF  Convective  heot  tronsfer  coefficient  (kW/(m2-C)) 


REAL  K.NSN 

Test  if  the  flow  is  laminar  or  turbulent  based  on  Reynolds  number. 

IF  (RE. GE. 2000.0)  THEN 

NSN=0.023»RE»*0.8*PR»*0.4 

ELSE 

NSN-3 . 66 
END  IF 

Colculote  the  convective  heot  transfer  coefficient  using  the 
Nusselt  number. 

HCOVF-NSN*K/D 

RETURN 

END 


SUBROUT  I NE  PRDPP (CARB . HYDR , OXYG . XNTR . SU LF . XA I R ) 


Combustion  products 
INPUTS: 

PT  Pressure  (otm) 
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XAIR 

Excess 

air  fraction 

CiRB 

Numbe  r 

of  moles  of  corbon 

H'tDR 

Numbe  r 

of  moles  of  hydrogen 

OXYG 

Numbe  r 

of  moles  of  oxygen 

XNTR 

Numbe  r 

of  moles  of  -nitrogen 

SULF 

Numbe  r 

of  moles  of  sulfur 

OUTPUTS: 

PC02 

Port i 0 1 

pressure  of  C02  (otm) 

PH20 

Partiol 

pressure  of  H20  (atm) 

PS02 

Po  r t i 0 1 

pressure  of  S02  (otm) 

PN2 

Partiol 

pressure  of  N2  (atm) 

P02 

Po  r t i a 1 

pressure  of  02  (atm) 

RATF 

Moss  rotio  of  air  to  fuel  (-) 

RPTF 

Moss  ratio  of  product  to  fuel 

(-) 

RWTF 

Moss  rotio  of  water  in  product 

t 0 

COMMON/PRODCT/PT  . PC02  . PH20 . PN2  , P02  . PS02  . RATF  , R.vTF . RPTF 

Moles  of  oxygen  and  nitrogen  to  be  odded 

R00=2  . •CA-’B+0 . 5*HYDR+2  . •SULF-OXYG 
RN0=3.76*R00 

Moles  of  C02 , H20.  02.  S02  ond  N2  in  products 

PC02=CARB 

PH2O=0.5»HYDR 

PO2=0.5»R00«XAIR 

PS02-SULF 

PN2»0.5»XNTR+0.5*RN0*(1 .+XAIR) 

Weight  of  fuel,  air,  and  products 

WT  F- 1 2 . • CARB+H YDR+ 16.*  0XYG4 14.*  XNTR+32 . • SU  LF 
WTA»(16. •R00+14. .RN0)*(1 .+XAIR) 

WTP=(44 . •PC02+1 8 . •PH20+32 . •P02+64 . •PS02+28 . .PN2) 
WTW*18. •PH20 


Compute  partiol  pressures  and  mass  ratios 

AMULT-PT/ ( PC02+PH20+P02+PS02+PN2 ) 

PC02-PC02*AMULT 

PH2CN=PH20*AKiRJLT 

PS02=PS02*AMULT 

P02-P02*AMULT 

PN2-PN2*AMULT 

RATF=WTA/WTF 

RPTF-WTP/WTF 

RWTF-WTW/WTF 

RETURN 

END 


SUBROUTINE  PRDPR(T,A»Ai(,AKK) 


Viscosity  ond  conductivity  of  combustion  product 


INPUTS: 

PC02 

Port i 0 1 

PH20 

Part iai 

PS02 

Port i 0 1 

PN2 

Partiol 

P02 

Port i 0 1 

OUTPUTS : 

M4A 

Product 

AKK 

Product 

pressure  of  C02  (otm) 
pressure  of  H20  (otm) 
pressure  of  S02  (otm) 
pressure  of  N2  (atm) 
pressure  of  02  (otm) 


dynomic  viscosity  (l<g/(m-s)) 
thermal  conductivity  (kW.(m^)) 
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COMMON/PRODCT/PT , PC02 . PH20 . PN2 . P02 . PS02 , RATF . RWTF . RPTF 
DIMENSION  X(1 1 ) ,A( 1 1 ) .B( 1 1 ) , AM(5) . AK(5) 

DATA  X/-1 0 . . 260 . , 440 . . 620 . . 800 . . 1 070 . . 
t1520. .2420. .3140. .3500. .4160./ 

DATA  A/0 . 0387 . 0 . 0553 . 0 . 0646 . 0 . 073 . 0 . 0806 ,0.0911. 

40.1062.0. 1320.0.1500.0.1583.0.1710/ 

DATA  B/0 .01287,0.01944.0. 02333 , 0 . 02692 . 0 . 03022 , 

40 . 03483 .0. 04178.0. 05348 .0.0612.0. 0646 . 0 . 0709/ 

DATA  AM/0 .889.0.685,1.123.0.964.0. 889/ 

DATA  AK/0 .925.0.906.1.037.0.983.0. 925/ 

NP=1  1 

AMM=PC02»AM( 1 )+PH20»AM(2)+P02«AM(3)+FN2»AM(4)+FSC2*AM(5) 
AKK=PC02 • AK ( 1 ) +PH2C • AK ( 2 ) +P02 • AK ( 3 ) +PN2 • AK ( 4 ) +PS02 • AK ( 5 ) 
TF=32.+1 .8*T 

AMM=(4. 1333E-04).AKW.TAB1 (NP.X.A.TF) 

AKK=( 1 .731E-03)»AKK»TAB1 (NP.X.B.TF) 

RETURN 

END 


FUNCTION  TAFF(HHV.TBS) 


Adiabatic  flame  temperature 


INPUTS: 

HHV  Higher  heating  value  of  fuel  (kJ/kg) 
RWTF  Weight  of  H20  per  kg  of  fuel  (kg/kg) 

RPTF  Weigth  of  product  per  kg  of  fuel  (kg/kg) 

PC02  Number  of  moles  of  C02 

PH20  Number  of  moles  of  H20 

P02  Number  of  moles  of  02 

PS02  Number  of  moles  of  S02 

PN2  Number  of  moles  of  N2 

TBS  Base  temperature  (C) 

OUTPUT : 

TAF  Adiabatic  flame  temperature  (C) 


COk*AON/PRODCT/PT , PC02 . PH20 , PN2 . P02 , PS02 , RAT  F . RWT  F , RPT  F 
ITR-0 

HW-HHV-2442 . •RWTF 
TAF  1 -0 . 995»HW/RPTF+TBS 
CP-CPF(TAF1 .TBS) 

TAF2-HW/RPTF/CP+TBS 

T-ABS(TAF1-TAF2) 

TAF1-TAF2 

ITR-ITR+1 

IF(ITR.LT.10  .AND.  T.GT.0.1)  THEN 
GOTO  100 
END  IF 
TAFF  -TAF1 
RETURN 
END 


FUNCTION  TABI(NP.X.Y.XI) 


Determine  Y(X1)  value  from  tabuloted  Y vs  X values 
INPUTS/OUTPUTS: 
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NP  Number  of  doto  points 

X Data  points  of  independent  variables 
Y Data  points  of  dependent  voriables 
TAB1  Output  dependent  varioble 

XI  Input  independent  variable 


DIMENSION  X(NP) .Y(NP) 

Set  out  the  range  Y1(X1)  values 

XMIN=X(1  ) 

XMAX=X(NP) 

IF(X1 . LE.XMIN)  THEN 
Y1=Y(1 ) 

ELSEIF(X1 .GT.XMAX)  THEN 
Y1=Y(NP) 

ELSE 

Interpolote  for  Y1(X1)  volue 


1 = 1 

IF  (XI .LE.X(I))  THEN 
DXT=X(I)-X(I-1) 
DYT=Y(I)-Y(I-1) 
R0XT=1 ./DXT 
DYODX=DYT*RDXT 
DX=X1-X( 1-1 ) 
Y1-Y(I-1  HDX*DY00X 
ELSE 
1 = 1 + 1 
GOTO  10 
END  IF 
END  IF 
TAB1-Y1 
RETURN 
END 
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APPENDIX  B Boiler  Simulation  Model  Setup 


al 


BLC22D  - Reduced  nuaber  of  reported  variables  to  BLC22C  880628 

SUPERBLOCK  1 
BLOCK  1 

UNIT  1 TYPE  62  - BOILER  MODEL  WITH  A TANKLESS  COIL 

UNIT  2 TYPE  63  - DOMESTIC  HOT  WATER  COIL 

UNIT  3 TYPE  64  - BOILER  BURNER  AND  WATER  PUMP  CONTROL 


UNIT  1 TYPE  62 

BOILER  MODEL  WITH  A TANKLESS  COIL 


1 


2 


3 


INPUTS : 


TEMPERATURE 

1 

- 

TSTK 

Stack  gas  temperature  (DE) 

TEMPERATURE 

2 

- 

TRA 

Room  air  temperture 

TEMPERATURE 

3 

- 

TOA 

Outdoor  air  temperature 

TEMPERATURE 

4 

- 

TBLW 

Boiler  water  temperature 

PRESSURE 

1 

- 

PSW 

Supply  water  pressure 

FLOW 

1 

- 

WSW 

Supply  water  flowrate 

TEMPERATURE 

5 

- 

TRW 

Return  water  temp. 

POWER 

1 

- 

QCW 

Heat  flowrate  to  domestic  hot  water 

CONTROL 

1 

- 

SWCH 

Control  switch  for  on/off 

OUTPUTS ; 

TEMPERATURE 

1 

- 

TSTK 

Stack  gas  temperature  (DE) 

TEMPERATURE 

4 

- 

TBLW 

Boiler  water  temp.  (DE) 

TEMPERATURE 

6 

- 

TSW 

Supply  water  temp. 

POWER 

2 

- 

QHWT 

Heat  gain  of  water 

POWER 

3 

- 

QSTK 

Heat  loss  from  stack 

POWER 

4 

- 

QINPUT 

Heat  input  rate 

POWER 

5 

- 

QLTNT 

Latent  heat  loss 

POWER 

6 

- 

QJAKT 

Heat  loss  thru  boiler  jacket 

POWER 

7 

- 

QSS 

Heat  flow  rate  to  boiler  water 

PARAMETERS : 

0.736000E-01  VGF 
0.175500  ASF 
0.000000  ARF 
0.381000  LGPF 
0.800000  EMISF 
0.776000E-01  VWB 
2 . 60000  AJAKT 
0.169290E-02  UJAKT 

1.00000  GARB 
1 . 84400  HYDR 

0 . 000000  OXYG 
0.000000  XNTR 
0.300000E-02  SULF 

45327.0  HFUEL 
0.839000  CPFUEL 
0.126600E-02  WFUEL 

26.5000  TFUEL 


Voltune  of  gas  in  the  boiler  fire-box  (m3) 
Boiler  fire-box  effective  radiation  heat  trans 
Boiler  fire-box  refractory  surface  area  (m2) 
Gas-path  length  in  the  boiler  fire-box  (m) 
Fire-box  surface  emissivity  value  in  fraction 
Volume  of  water  in  the  boiler  (m3) 

Boiler  jacket  surface  area  (m2) 

Boiler  jacket  U-factor  (kW/m2-C) 

Atomic  ratio  of  carbon  in  fuel 
Atomic  ratio  of  hydrogen  in  fuel 
Atomic  ratio  of  oxygen  in  fuel 
Atomic  ratio  of  nitrogen  in  fuel 
Atomic  ratio  of  sulfur  in  fuel 
Fuel  higher  heating  value  (kJ/kg) 

Fuel  specific  heat  value  (kJ/kg-C) 

Fuel  supply  rate  (kg/s) 

Fuel  temperature  (C) 
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0.200000 

XAIR 

0.790000 

EFFYSS 

95.0000 

TBLWS 

320.000 

TSGS 

0.925000 

ICFGON 

3.14000 

ICFGOf 

1.38000 

ICFW 

3.00000 

DTFG 

0.231250 

ICFGNS 

0.462500 

ICFGFS 

30.0000 

DTBW 

Execess  combustion  air  (-) 

Boiler  full  load  efficiency  (-) 

Boiler  water  temp,  at  steady  state  (C) 

Stack  gas  temp,  at  full  load  (C) 

Integration  control  for  on-period  stack  gas 
Integration  control  for  off-period  stack  gas 
Integration  control  for  boiler  water  temperatu 
Short  start  period  for  burner  on/off  cycle  (s) 
Integration  control  for  short  start  on-period 
Integration  control  for  short  start  off-period 
Time  interval  for  updating  QSS  (s) 


UNIT  2 TYPE  63 

DOMESTIC  HOT  WATER  COIL 


1 

INPUTS: 

TEMPERATURE 

4 - 

TBLW 

; Boiler  water  temperature 

FLOW 

2 - 

WCW 

; Water  flowrate  thru  the  coil 

TEMPERATURE 

7 - 

TICW 

; Coll  inlet  water  temp. 

2 

OUTPUTS: 

TEMPERATURE 

8 - 

TOCW 

; Coil  outlet  water  temperature 

POWER 

1 - 

QCW 

; Heat  flowrate  to  the  coil  water 

3 

PARAMETERS: 

0.193000E- 

■01  DCOIL 

: Diameter  of  the  coil  (m) 

4.20600 

LCOIL 

: Length  of  the  coil  (m) 

UNIT 

3 TYPE 

64 

BOILER  BURNER  AND 

WATER  PUMP  CONTROL 

1 

INPUTS: 

TEMPERATURE 

10  - 

TBURN 

: Burner  control  temperature 

TEMPERATURE 

4 - 

TBLW 

: Boiler  water  temperatxire 

TEMPERATURE 

9 - 

TLOAD 

: Load  temperature 

CONTROL 

2 - 

CTHERM 

: Thermostat  control  indicator  1/0  : 

FLOW 

2 - 

WCW 

: Water  flowrate  thru  the  coil 

2 

OUTPUTS: 

CONTROL 

1 - 

SWCH 

: Burner  control  on/off  signal 

FLOW 

1 - 

WSW 

: Boiler  circulating  water  flowrate 

3 


(- 


PARAMETERS: 

85.0000 

TBURON 

98.9000 

TBUROF 

70.0000 

TPMPON 

80.0000 

TPMPOF 

0.286100 

WPUMP 

O.lOOOOOE-03  RATEOF 


Boiler  bumer-on  temp.  (C) 

Boiler  bumer-off  temp.  (C) 

Boiler  water  pump-on  temperature  (C) 

Boiler  water  pump-off  temperature  (C) 
Circulating  water  flowrate  (kg/s) 

Minimum  water  flow  rate  ratio  when  pump  off 


Initial  Variable  Values: 
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PRESSURE 

1 -> 

0.000000 

(kPa) 

FLOW 

1 -> 

0.000000 

(kg/s) 

FLOW 

2 -> 

0.000000 

(kg/s) 

TEMPERATURE 

1 -> 

315.000 

(C) 

TEMPERATURE 

2 -> 

26.5000 

(C) 

TEMPERATURE 

3 -> 

17.0000 

(C) 

TEMPERATURE 

4 -> 

30.0000 

(C) 

TEMPERATURE 

5 -> 

57.3800 

(C) 

TEMPERATURE 

6 -> 

30.0000 

(C) 

TEMPERATURE 

7 -> 

14.0000 

(C) 

TEMPERATURE 

8 -> 

28.0000 

(C) 

TEMPERATURE 

9 -> 

25.0000 

(C) 

TEMPERATURE 

10  -> 

80.0000 

(C) 

CONTROL 

1 -> 

1.00000 

(-) 

CONTROL 

2 -> 

0.000000 

(-) 

POWER 

1 -> 

0.000000 

(kW) 

POWER 

2 -> 

0.000000 

(kW) 

POWER 

3 -> 

0.000000 

(kW) 

POWER 

4 -> 

0.000000 

(kW) 

POWER 

5 -> 

0.000000 

(kW) 

POWER 

6 -> 

0.000000 

(kW) 

POWER 

7 -> 

0.000000 

(kW) 

simulation  Error  Tolerances: 

1 RTOLX=  O.lOOOOOE-03  ATOLX=  O.lOOOOOE-04 
XTOI^  0.200000E-03  TTIME=  1.00000 

SUPERBLOCK  1 

2 FREEZE  OPTION  0 SCAN  OPTION  0 


The  following  are  Boundary  Variables  in  the  simulation: 


FLOW  2 
TEMPERATURE  9 
TEMPERATURE  10 
CONTROL  2 


The  following  are  the  reported  variables: 


SUPERBLOCK  1 REPORTING  INTERVAL  120.000 

TEMPERATURE  1 

TEMPERATURE  4 

TEMPERATURE  8 

CONTROL  1 

POWER  1 

POWER  2 

POWER  4 
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APPENDIX  C Boundary  Data  File  Generation  Program  - CRBND 


CRBND3 

: Gene  rot i 

on  of  boundary  doto 

f i 

1 e for 

the 

bo  i 

1 1 e r mode  1 

when  domestic  hot  woter  is 

dr 

Own  with 

epoce 

heo  t i ng . 

This  prog 

rom  generotes  o set 

of 

doto  f i 

les 

or 

0 file 

for  spoce 

ioods  ot  0 given  dome 

Stic  hot 

water 

d row . 

August 

5.  1988 

Cheol  Pork 

PROGRAM  CRBND3 


PARAMETER  (MAXKt«15) 

REAL  XPUMP(MAXN) .CYCLERT(MAXN) .NCYCLE 

INTEGER  TSTOP,TPS.TPON(MAXN) ,TPOFF(MAXN) ,NCYFWP(MAXN) 

CHARACTER  BNDFILE*12.  FILEXT(MAXN>3,STEPPMP*1 1 . ERRFl lE*12 

character  ST epout* 12, answer *1 


DATA 

t 

DATA 

t 


XPUMP/e  02.  0.05,  0.10.  0.15.  0.20.  0.225.0.30. 

0.40.  0,50.  0.60,  0.70,  0.80.  0.90.  0.95.  0 975/ 
CYCLERT/0.4,  0.6,  1.20.  1.40.  1.70.  1.80,  2.20. 
2.45.  2.5.  2.35.  2.20.  1.70.  0.90.  0.50.  0.40/ 


DATA  TPS/ 12/ 


DATA 

* 


FILEXT/’002’ . •005’ . *010’ . •015’ . ’020' . •023’ . •030' . 
•040' . •050’ . •060^ . ^070^ . •080^ . •090^. •095^ . •098^/ 


PRINT  Is  this  run  STEP  1 or  STEP  2 ?• 
READ  ‘.ISTEP 

PRINT  •. • ISTEP  - • . ISTEP 


Dstermins  the  pump  on/off— times  ond  the  number  of  cycles  per 
the  period. 

PRINT  1000 
T STOP-5000+TPS*  3600 
DO  10  1-1 ,MAXN 
NCYCLE-CYCLERT ( I ) .TPS 
TPON(  I )-3600.XPULP(  I )/CYCLERTn  ) 

TPOFF( I )-3600.TPS/NCYCLE-TPON( I ) 

IF(NCYCLE-INT(NCYCLE)  .GT.  0.0)  THEN 
NCYP1mP(  I )-INT(NCYCLE+1 .0) 

ELSE 

NCYPWP( I )-NCYCLE 
END  IF 

PRINT  2000.  I.XPUIyP(I).CYCLERT(I).NCYPMP(I).TPON(I).TPOFF(I). 
k TSTOP 

10  CONTINUE 

C Creote  o boundory  doto  files  with  different  file  extension 

IFdSTEP  .EQ.  1)  then 
CALL  STEPINA 
ENDIF 

PRINT  .,•  Are  you  processing  o tingle  file  ? (N) • 

READ  (..FMT-^CAd^)  ANSWER 

IF(ANSWER  .EO.  'Y^  .OR.  ANSWER  .EQ.  ‘y’)  THEN 
PRINT  .,•  Whot  is  the  index  number?^ 

READ  «.  INDEX 
I-INDEX 

STEPPMP-^STEPPW^.  •//FILEXTHI 
STEPOUT--STEPOUT.  •//FILEXTH) 

BNDF I LE- • BNDE000 . • //F I LEXT  f I ) 

ERRF1lE-^ERRE000. •//FILEXT(l) 


72 


ooooo  oo 


1F(ISTEP  .EQ.  1)  THEN 

CALL  STEPIN(TPON( I ) ,TPOFF( I ) , NCYPMP( 1 ) . STEPPMP) 

CALL  STEPCOM(STEPOUT) 

PRINT  •,*  Edit  STEPOUT; • ,FI LEXT(I ) , • ond  perform  STEP  2 ' 
ELSE1F(ISTEP  .E0.2)  THEN 

CALL  STEPBND(BNDFI LE . STEPOUT . NUMCOL) 

CALL  CHEKBND(BNDFI LE , ERRFI LE . NUMCOL) 

PRINT  Edit  BNDE000. ' .FILEXT(I) 

PRINT  based  on  ERRE000 . ’ . FI LEXT( I ) 

PRINT  ■■'nil 

PRINT  • 

END  IF 
ELSE 

PRINT  •,*  NOTE;  Excluded  ore  Index  1 ond  15.  ’ 

DO  20  I-2.MAXN-1 
C«  DO  20  1=1 .MAXN 

STEPPMP=’ STEPPMP. 

STEPOUT»’ STEPOUT.  '//ZIIZXTU) 

BNDFILE=’BNDE000.  '/AILEXTO) 

ERRFILE-'ERRE000. •//FltZXl(l) 

IFdSTEP  .EO.  1)  THEN 

CALL  STEPIN(TPON( I ) .TPOrF( I ) .NCYPMP( I ) . STEPPWP) 

CALL  ST EPCOM( STEPOUT) 

PRINT  *.•  Edit  STEPOUT. ' .FI LEXT( I). * ond  perform  STEP  2 ‘ 
ELSEIF(ISTEP  .EQ.2)  THEN 
CALL  STEPBND(BNDFI LE . STEPOUT . NUMCOL) 

CALL  CHEKBND(BNDFI LE . ERRFI LE , NUMCOL) 

PRINT  Edit  BNOE000. • ,FILEXT(I) 

PRINT  •.*  based  on  ERRE000. ’ .FILEXT(I) 

PRINT  ■ I 

PRINT  • 

END  IF 

20  CONTINUE 

ENDIF 

1000  F0RMAT(/79C*’)/T3.  T .T8. ‘Xpump’ .T18. ’Cycle  Rote’.T30. 

k • # cycles’. T42.’Tp. on’. T52.’Tp. off ’.T62.’T8top’/79(’.’)//) 
2000  FORMAT(I5.2F10.3.4I10) 

STOP 

END 


SUBROUTINE  STEPINA 

STEPINA  : Generotion  of  input  doto  to  STEPBNO.FTN 
when  domestic  hot  woter  is  drown. 


IMPLICIT 
CHARACTER* 1 
REAL 


INTEGER(A-Z) 

YESNO 

WCWON.WCWOFF.US 


DATA  TBURN/100/.CTHERM/1/ 

DATA  YESNO/’y’/.TLDON.TLDOFF/68.85/ 


OPEN ( 8 . F I LE- ’ ST  EP I NA2 . DAT ’ ) 


PRINT  *.’  Enter  the  totol  omount  of  woter  drown  in  Gol’ 
READ  *.  US 


PRINT  *.’  How  ffiony  cycles  in  6-hour  period?’ 
READ  *,  NCYCLE 


PRINT  •.’  Stort  time  of  drow?’ 
READ  *.  TSTART 
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WCWON=0. 1892 
WCWOFF-=0.0 

DURATN=US*3 . 7854/(WCW0N*NCYCLE) 

PRINT  DURATION  OF  DRAW  -=’.DURATN 

PRINT  DRAW  RATE  (KG/S)  -’.WCWON 

DO  10  I=1.NCYCLE 

T I MEON-TST ART+ ( I - 1 ) • 3600*  6/NCYC  LE 
T I MEOF«=T  I MEON+DURATN 


WRITE(8.1000) 
WRITE(8.1100) 
WRITE(8. 1 100) 
WRITE(8.1000) 
WRITE(8.1100) 
WRITE(8,1100) 


YESNO 

T IMEON , WCWOFF . TLDOFF , TBURN . CTHERM 
TIMEON.WCWON,  TLDOFF . TBURN . CTHERM 
YESNO 

TIMEOF.WCWON,  TLDOFF, TBURN. CTHERM 
TIMEOF, WCWOFF. TLDOFF. TBURN. CTHERM 


10  CONTINUE 


YESNO’ n' 


WRITE(8. 1000)  YESNO 

1000  FORMAT (A1) 

1100  FCRMAT(I10.F10.4.3I5) 


RETURN 

END 


SUBROUT I NE  ST EP I N ( TPON . TPOFF . NCYPMP . ST EPPMP ) 

STEPIN  : Generotion  of  input  doto  to  STEPBND.FTN 
when  the  woter  punp  is  cycled. 


IMPLICIT  INTEGER(A-Z) 
CHARACTER  STEPPlyP*  1 1 
CHARACTER* 1 YESNO 


C* 


DATA  TSTART/5000/ . WCW/0/ . TBURN/1 00/ . CTHERM/1/ 

DATA  TLDON.TLDOFF/68.85/ 

NAMELIST  /NAM2/  TPON. TPOFF. NCYPMP 


OPEN(7.FILE-STEPPMP) 


C«  PRINT  NAM2 

OO  10  1-1 .NCYPMP 
YESNO’ y’ 

TIMEON-TSTART+( 1-1 )• (TPON+TPOFF) 
TIMEOF-TSTART+I •TPON+( 1-1 ) •TPOFF 


WR1TE(7.1000) 

WRITE(7.1100) 

WRITE(7.1100) 

WRITE(7.1000) 

WRITE(7.1100) 

WRITE(7.1100) 


YESNO 

T I MEON . WCW , T LDOFF . TBURN .CTHERM 
TIMEON.WCW.TLDON  .TBURN. CTHERM 
YESNO 

TIMEOF. WCW. TLDON  .TBURN. CTHERM 
T I MEOF . WCW . T LDOFF . TBURN . CTHERM 


10  CONTINUE 


YESNO’ n’ 

TSTOP-T I MEOF+TPOFF 


WRITE(7,1000)  YESNO 
WRITE(7,1100)  TSTOP 

1000  FORMAT(AI) 

1100  FORMAT(I 10.415) 
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RETURN 

END 


SUBROUTINE  STEPCOM(STEPOUT)' 

STEPCOM  : Combining  outputs  of  STEPIN.FTN  for  space  heating 
and  domestic  hot  water  drew 


PARAMETER  (MAXDAT-500) 

CHARACTER* 1 YESNO 
CHARACTER *80  HEAD 

CHARACTER  STEPPMP*! 1 .STEP0UT*1 1 .FILEXT*3 
REAL  WCWS(MAXDAT,3) .WCWD(MAXDAT,3) 

INTEGER  TIMES(MAXDAT) ,TLDS(MAXDAT.3) .TLDD(MAXDAT.3) 

INTEGER  TBURN.CTHERM.TIMED(MAXDAT) ,TIME(MAXDAT) . 

IWCWS(MAXDAT.3) .TSTOP 

0PEN(9,FILE=STEP0UT) 

OPE--'(10,FILE-’STEPLEAD.DAT’  ) 

REV»  NO  7 
REWIND  8 
REWIND  10 


Read  File  10  and  write  on  File  9 ( In i t i a I i zot i on  heading  ) 


L-1 

5 READ(10. 1100, END-7)  HEAD 

WRITE{9.1100)  HEAD 

GOTO  5 


• Read  File  7 ( Output  for  space  heoting) 

7 1-1 

10  REAO(7.1000)  YESNO 

IF(YESNO  .EQ.  ’y’)  THEN 

READ(7.*)  TIMESn).IWCWS(I.1).TLDSfI.1).TBURN.CTHERM 
READ(7.*)  TIMES(I).IWCWS(I.2).TLDS(I.2).TBURN.CTHERM 
I-I+1 
GOTO  10 
ENDIF 

READ (7.*)  TSTOP 
NUM0AT7-I-1 

PRINT  *,'  Totol  number  of  time  chonges  of  file  7 - NUM0AT7 

• Reod  File  6 ( Output  for  domestic  hot  woter  draw) 

J-1 

20  READ(8.1000)  YESNO 

1F(YESN0  .EG.  'y')  THEN 

READ(8.*)  TIMED(J).WCWDfJ.1).TLDD(J.1) .TBURN.CTHERM 
READ(8.*)  T1MED(J),WCWD(J.2).TLDD(J.2).TBURN.CTHERM 
J-J+1 
GOTO  20 
ENDIF 

NUMDAT8-J-1 

PRINT  Total  number  of  time  changes  of  file  8 - NUMDAT8 

• Combining  two  file  with  o kind  of  sorting 

TIMEMK-TIMED(I) 

1 START-1 

JJ-1 

KK-0 

YESNO-'y* 

25  DO  40  1-ISTART.NUMDAT7 

1F(TIMES(I)  .LT.  TIMEkK)  THEN 
K1-KK+I-ISTART+1 
T1ME(K1)-TIMES(I) 
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wcws(i.i)-iwcws(i,i) 

WCWSd  .2)=IWCWS(1  .2) 

WRITE(9, 1000)  YESNO 

WRITE(9. 1200)  TIMES(I) .WCWSn , 1 ) . T LDS( I . 1 ) . TBURN . CTHERM 
WRITE(9. 1200)  TIMES(l) .WCWSO .2) .TLDS(I . 2 ). TBURN , CTHERM 
C«  PRINT  .,'1.  K1.  TIME(K1)-M  ,K1  .TIME(KI) 

ELSE 

DO  30  J=1 .NUMOAT8 
IF(TIMED(J)  .GT.  TIMES(I-I)  .AND. 
k TIMED(J)  .LE.  TIMES(I))  THEN 

K2-K1+J-JJ+1 
TIME(K2)-TIMED(J) 

WRITE(9. 1000)  YESNO 

WR1TE(9. 1200)  TIMED(J) .WCWD(J , 1 ) .TLDD(J . 1 ) .TBURN. CTHERM 
WRITE(9.1200)  TIMED(J) .WCWO(J. 2) .TLDD(J.2) .TBURN. CTHERM 
C*  PRINT  •.•J.  K2.  TIME(K2)-‘ .J.K2.TIME(K2) 

ELSEIF(TIMED(J)  .GT.  TIMES(I))  THEN 

jj-j 

TIMEMK-TIMED(J) 

GOTO  50 
END  IF 


30 

CONTINUE 

END  IF 

40 

CONTINUE 

50 

KK=K2 

ISTART-I 

I F(  I START  .LT.  NUM-.AT7)  THEN 

GOTO  25 
END  IF 


DO  55  I-1.NUMOAT7 
IF(TIMES(I)  .GE.  TIMEMK)  THEN 
K3-NUMDAT8+I 
TIME(K3)-TIMES(I) 

WCWS(I .l)-IWCWS(I.I) 

WCWS(I .2)-IWCWS(I .2) 

WRITE(9.1000)  YESNO 

WRITE(9. 1200)  TIMES(I) .WCWS(I . 1 ) .TLDSf I , 1 ) .TBURN. CTHERM 
WRITE(9.1200)  TIMES(I).WC1WS(I .2).TLDS(I.2).TBURN.CTHERM 
C*  PRINT  K3.  TIME(K3)  . I .K3.TIME(K3) 

END  IF 

55  CONTINUE 

YESNO-' n’ 

WRITE(9.1000)  YESNO 

• Write  output 

C*  DO  60  1-1 .NUMDAT7+NUMOAT8 

C.  PRINT  I.  TIME(I)  I.  TIME(I) 

C*  WRITE(9.1100)  I.  TIMEO) 

C»60  CONTINUE 

CLOSE(9) 

1000  FORMAT (A1) 

1100  FORMAT (A80) 

1200  FORMAT(I10.F10.4.3110) 

RETURN 

END 


SUBROUT  I NE  ST  EPBND ( BNDF I LE . STEPOUT . NUMCOL ) 
STEPBND  : Stepwise  changes  in  Boundary  variables 


t0  t1  t2  t3  t4  t5  t6  t7  t8  t9 

t 

I 

The  step  change  occurs  at  t1. 
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NUMSUB  ; Number  of  subintervols  in  on  interval  (=5) 

NUWINT  ; Number  of  intervols  chosen  (=9) 

DELTSUB;  Sire  of  eoch  subinterval  in  seconds 
DELT  : Sire  of  eoch  interval  in  seconds 
NUMCOL  : Number  of  columns  of  Y— voriobles 

In  the  figure  shown  above,  there  ore  9 intervols,  of  which 
sire  is  DELT  seconds,  DELT  - NUMSUB*DELTSB . 


PARAMETER 

CHARACTER 

REAL 

REAL 

REAL 

INTEGER 


(MAXCOL-30 . MAX  I NT-20 , MAXSUB-1 0 ) 

ANSWER* 1 , BNDF I LE* 1 2 , FI LEXT.3 , STEPOUT* 1 2 
Y1  (MAXC0LKY2(MAXC0L)  ,DELY(MAXCOL)  ,YY(MAXCOL) 
TIME(0:MAXINT) ,TIME1 (0.MAXINT,0:MAXSUB) 
Y(0;MAXINT,MAXCOL) 

TIME2 


OPEN(9,  FILE-STEPOUT) 
OPEN(11 ,FILE=BNDFILE) 
CLOSEdI  ,STATUS-’DELETE') 

OPEN(11 ,FILE»BNDFILE) 
REWIND  9 


PRINT  *,  ' A boundory  doto  is  being  generated  ’ 

C*  PRINT  *,  ’ What  is  the  number  of  columns  for  Y— volues  ?' 

READ(9,*)  NUMCOL 

C*  PRINT  *,'  Enter  Subinterval  sire  in  seconds’ 

C*  PRINT  •,*  Number  of  Intervols?,  ond  Number  of  Subintervals' 

READ(9,*)  DELTSB,NLA4I NT. NUMSUB 
IF(NUMSUB  .LT.  4)  THEN 

PRINT  *.  • NUMBER  OF  SUBINTERVALS  MUST  BE  MORE  THAN  4 ■’ 

STOP  • MOOIFIY  YOUR  INPUT  DATA  AND  TRY  AGAIN  ’ 

END  IF 

C*10  PRINT  •,*  Enter  Step  chonge  time,  ond  Y-volues  before  change’ 
10  CONTINUE 

READ(9,*)  TSTEP.(Y1( I). 1-1 .NUMCOL) 

C*  PRINT  Enter  Step  change  time,  ond  Y-volues  after  chonge’ 

READ(9.*)  TSTEP.(Y2( I). 1-1 .NUMCOL) 

OELT-NUMSUB*DELTSB 

TIME(0)-TSTEP-OELT 

DO  20  J-0,NUMINT 

TIME(J)-TIME(0)+REAL(J)*DELT 
DO  20  K-0. NUMSUB 

TIME1(J.K)-TIME(J)+REAL(K)*DELTSB 
20  CONTINUE 

DO  30  1-1 .NUMCOL 

DELY( I )-(Y2( I )-Y1 ( I ) )/(NUMINT-1 ) 

Y(0.I)-Y1(I) 

DO  30  J-0.NUMINT-1 

Y(J.I)-Y(0.I)+REAL(J)*OELY(I) 

30  CONTINUE 

DO  60  J-0.NUMINT-1 
DO  40  1-1 .NUMCOL 
40  YY(I)-Y(J.1)+1 .0E-10 

DO  50  K-0. NUMSUB 

TIME2»TIME1(J.K)+1 .0E-10 

WRITE(11 ,1000)  TIME2.(YY(L).L-1 .NUMCOL) 

50  CONTINUE 

60  CONTINUE 

C*  PRINT  •.‘Continue?’ 

READ  (9.FMT-’(A1)’)  ANSWER 

1F(ANSWER  .EO.’Y’.OR.  ANSWER  .EQ.  ’y’)  GOTO  10 
1000  FORMAT(1X.I7.F10.4.6F10.2) 
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RETURN 

END 


SUBROUT  1 NE  CHEKBND ( BNDF I LE . ERR F I LE . NUMCO L ) 

CHEKBND  : Check  the  boundory  date  file  for  proper  sequence 
of  the  time  dota 


PARAMETER  (MAXCOL-30) 

CHARACTER  ERRFI LE» 12 . BNDFI LE» 1 2 
REAL  YY(MAXCOL) 

INTEGER  TIMENEW.TIMEOLD.TIME2 

OPEN(12,FILE-ERRFILE) 

CLOSE(12.STATUS-*DELETE‘) 

0PEN(12.FILE-ERRFILE) 

OPEN(11 .FILE=BNDFILE) 

REWIND  11 

TIMEOLD-0 

10  READ(11 .10e0,END=20)  TIME2. (YY(L) , L“1 .NUMCOL) 

TIMENEW-TIME2 

1F(TIMENEW  .LT.  TIMEOLD)  THEN 

C*  PRINT  Time  is  less  thon  the  previous  volue  ’ 

PRINT  *.•  ERROR:  ’ .TIMENEW. TIMEOLD 
WRITE(12.2000)  TIMENEW. TIMEOLD 
ELSE 

TIMEOLD-TIMENEW 
END  IF 
GOTO  10 

20  CONTINUE 

1000  FORMAT(1X.I7.F10.4.6F10.2) 

2000  FORMAT(2X.2I12) 

RETURN 

END 
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APPENDIX  D Integration  Routine  for  Energy  Values  - ENERGY 


ENERGY  : Integrotion  by  using  traperzoidol  method 

Step  change  is  considered  for  the  SB*. DAT  files 

Moy  10.  1988  C.  P. 

Revised  : July  15,  1988 


NCOL  : Number  of  doto  columns  per  o record 

NROW  : Number  of  data  rows  in  o whole  file 

ISEL  : Index  number  of  couimn  of  selected  data 

NRPR  : Number  of  rows  per  o record 

TBEGIN  : Begining  time  of  integration  (s) 

TEND  : Ending  time  of  integrotion  (s) 


PROGRAM  ENERGY 

PARAMETER  (MAXCOL-30.MAXROW- 15000) 

REAL  X(MAXCOL).Y(MAXROW),T(M  <ROW) 

CHARACTER  INPFIL*20.XX*80 

C DATA  NRPR/4/ 

DATA  DTM IN/0.1/ 

C*  PRINT  •,*  Doto  Fi I e Nome  ?’ 

C*  READ(*.FMT-‘(A20)*)  INPFIL 

C*  PRINT  1000,  INPFIL 

C*1000  FORMAT(10C*').A20.5X.10C*’)) 

C*  0PEN(7,FILE-INPF1L) 

5 0PENf8,FlLE-’ENERGY8.0UT’) 

READ(8.FMT-'(A80)* .END-7)  XX 
GOTO  5 

7 PRINT  •.*  Totol  number  of  columns  ond  rows  7* 

READ  •.  NCOL.  NROW 

PRINT  *,  ’ Number  of  rows  per  o record  in  the  fife  ?’ 
READ  • . MtPR 

VIIRITE(B.1000) 

10  REWIND  7 

PRINT  •,*  Selected  channel  number  ?’ 

READ(*.  •. END-50)  ISEL 

PRINT  e,  * Selected  chonnel  number  ISEL 

PRINT  •.*  Begining  time.  Ending  time  (s)  ?' 

READ  •.  TBEGIN,  TEND 

PRINT  TBEGIN  - • .TBEGIN,  * TEND-  ' ,TEt© 

PRINT  •,*  Maximum  time  intervol  ond  change  in  volue  * 
READ  •.  DTREF,  DYREF 

PRINT  •.*  DTREF  -’.DTREF, • DYREF  .DYREF 
I COUNT-1 

20  READ(7.w. END-30)  (X( J) . J-1 .NCOL) 

TIME  - X(1) 

1F(TIME  .GE.  TBEGIN  .AND.  TIME  .LE.  TEND)  THEN 
T(ICOUNT)-X(1) 

Y(ICOONT)-X(ISEL) 

ICOUNT-ICOUNT+I 
ELSEIF(TIME  .GT.  TEND)  THEN 
GOTO  30 
ENDIF 

1F(IC0UNT-1  .LT,  NROW/NRPR)  THEN 
GOTO  20 
ENDIF 

30  NCOUNT-ICOUNT-1 

IF(^COONT  .GE.  MAXROW)  THEN 

STOP  • MAXROW  IS  GREATER  THAN  15000  ' 
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END  IF 


SUK*-0 . 0 

DO  40  I-1,NCOUNT 

C Troperoidol  integration 

IF(I  .GT.  1)  THEN 

C Step  chonge 

IF((T(I)-T(I-1))  .GE.  DTREF  .AND. 
k (Y(I)-Y(I-I))  .GE.  DYREF)  THEN 

TR-T(I)-DTMIN 
S-Y(I-1).(TR-T(1-1)) 

ELSEIF((T(I)-T(I-1))  GE.  DTREF  .AND. 
ft  (Y(l-l)-Y(I))  .GE.  DYREF)  THEN 

TR-T(1-1)+DTMIN 
S-Y(I).(T(I)-TR) 

ELSE 

S=(Y(I)+Y(I-1)).(T(I)-T(I-1))/2.0 
END  IF 
ELSE 
S=*0.e 
END  IF 

C»  PRINT  I - •.  1.*  S -'.S 

SUM«SUVt+S 
40  CONTINUE 

AVE-SUM/ ( T END-TB  EG I N ) 

PRINT  TOTAL  - ’.SUM 

PRINT  AVERAGE-  AVE 

PRINT  • 

PRINT  

WRITE(8.2000)  ISEL.TBEGIN, TEND. SUM. AVE. DTREF. DYREF 

IF(ISEL  .EO.  5)  THEN 
TBON-SUM 

TBOFF-TEND-TBEG I N-TBON 
ELSEIFdSEL  .EQ.  7)  THEN 
OHWT-SUM 

ELSEIFdSEL  .EQ.  8)  THEN 
OINP-SUM 

ELSEIFdSEL  .EQ.  6)  THEN 
OCW-SUM 
END  IF 

GOTO  10 

50  XBURN-  TBON/(TBON+TBOFF) 

EF-  Qa«/QINP«100. 

PRINT*.'  TBON  -• .TBON. ' TBOFF  .TBOFF 
PRINT  *.'  XBURN  .XBURN.  ' EF-'.EF 

EFFU-  QHWT/QINP*100. 

PRINT  *.'  XBURN  -• .XBURN.  * EFFU  -’.EFFU 

1000  FORMAT(79('-')) 

2000  FORMATd5.2F7.0.F12.1.3Fl0.4) 

CLOSE(8) 

STOP  • END  OF  JOB  ' 

END 
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APPENDIX  E B-Sline  routine  - SPLINE 




c 

C SPL  : Interpolate  the  doto  points  using  B-spline  method 
C 

C 

c 

C August  22,  1968  Cheol  Pork 
C 

c 

c 

PROGRAM  SPL 

PARAMETER  (NDATA-100.NSETMAX-10) 

CHARACTER  TITLE(2)*80 

DIMENSION  X(NDATA  .Y(NDATA) .FDP(NDATA) .XX(NDATA) .YY(NDATA) 
DIMENSION  YNEW(NOATA.NSETMAX) 

DATA  N/12/.NSET/6/ 

OPEN(7,FILE-’SPL.INP') 

OPEN(8.FILE--SPL.X‘) 

OPEN(9.FILE-’SPL.OUT’) 

CLOSE(9. STATUS- ’DELETE') 

OPEN(9.FILE-'SPL.OUT') 

REWIND  7 
REWIND  8 
REWIND  9 


C Reod  the  reference  data  file. 

DO  60  KK-1.NSET 

READ (7. 2000. END-999)  TITLE(I) 

READ(7.2000)  TITLE(2) 

DO  10  1-1 .N 

READ(7.*)  IXPUMP.XBURNER,  X(1).Y(1) 
ie  CONTINUE 


C Calculate  the  derivotives  ot  eoch  of  point  (x,y) 

. CALL  SPLINE(N.X.Y.FOP) 

C Reod  the  input  file  of  XX. 

DO  30  J-1  .NDATA 
REA0(8.«.END-4e)  XX(J) 

30  CONTINUE 
40  NN^J-1 

REWIND  B 

DO  50  K-1.NN 

C Interpolote  for  given  date. 

CALL  SPEVAL{N.X.Y.FDP.XX(K).YY(K)) 

PRINT  •.  K.  XX(K).  YY(K) 

YNEW(K.KK)-YY(K) 

50  CONTINUE 
60  CONTINUE 

DO  70  K-1  .NN 

WRITEt9.1000)  XX(K).(YNEW(K,KK).KK-1 .NSET) 

70  CONTINUE 

1000  FORMAT (BF1 0.2) 

2000  FORMAT (AB0) 

699  STOP 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 


ie 


c 

c 

c 

c 


20 


30 

c 

c 

c 

c 

c 


END 


SUBROUTINE  SPLINE(N , X . Y . FDP) 


SPLINE  ; computes  the  second  derivotives  needed  in  cubic 

spline  interpolation.  The  original  program  was  written  by 
J.  M.  Ferriger  Ref.[l].  A little  modification  is  made. 

July  18.  1984  C.P. 


N: 

X: 

Y: 


FDP; 


Number  of  data  points 

Array  containing  the  values  of  the  independent  variable 
(Assume  to  be  in  oscending  order) 

Arroy  contoining  the  volues  of  the  function  at  the  data 
points  given  in  the  X orray 

Output  orray  which  contoins  the  second  derivotives  of 
the  interpolating  cubic  spline. 


REFERENCE: 

[l]  Joel  H.  Ferziger,  "Numericol  Methods  for  Engineering 

App I i cot i on . " John  Wiley  4 Sons,  1981,  pp. 17-18. 


REAL  LAMDA 

parameter  (NMAX-100) 

DIMENSION  X(NMAX) .Y(NMAX) .A(NMAX) .B(NMAX) ,C(NMAX) .R(NMAX) . 
4 FDP(NMAX) 

Compute  the  coefficients  ond  the  RHS  of  the  equotions. 

This  routine  uses  the  conti lever  condition.  The  parameter 
LAVOA  is  set  to  1.  But  this  con  be  user-modified. 

A,B,C  are  the  three  diagonols  of  the  tridiogonol  system, 
and  R is  the  right  hond  side. 

LAMDA  - 1 . 

C(1)-X(2)-X(1) 

DO  10  1-2. N-1 

cn)-x(i+n-x(i) 

An)-c(i-i) 

B(I  -2..(A(I)+C(1)) 

R(I)-6..((Y(1+1)-Y(I))/C(I)-(Y(I)-Y(I-1))/C(I-1)) 

CONTINUE 

B(2)-8(2)+LAMDA.C(1) 

B(N-1 )-B(N-1)+LAM0A*C(N-1) 

Tridiogonol  solver  subroutine. 

But  the  notoion  is  clumsy  so  we  will  solve  directly. 


DO  20  1-3, N-1 
T-A(I)/B(I-1) 

B(I)-B(I)-T.C(I-1) 

R(I)-R(I)-T*R(I-1) 

CONTINUE 

FDP(N-1 )-R(N-1 )/B(N-1 ) 

DO  30  1-2, N-2 

FDP(N-I)-(R(N-I)-C(N-I).FDP(N-I+1))/B(N-I) 

CONTINUE 

FDP(1 )-LAMDA.FDP(2) 

FDP(N)-LAMDA.FDP(N-1 ) 

RETURN 

END 


SUBROUTINE  SPEVAL(N.X.Y.FDP.XX.F) 
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speval 


evoluotes  the  cubic  spline  for  given  the  derivotives 
computed  by  subroutine  SPLINE. 


XX: 

F; 


Value  of  independent  vorioble  for  which  on  interpoloted 

value  is  reqested 

The  interpolated  result 


PARAMETER  (NMAX=100) 

DIMENSION  X(NMAX) ,Y(NMAX) ,FDP(NMAX) 

Find  the  proper  intervol. 

DO  10  1-1 .N-1 
IF(XX.LE.X(I+1))  GOTO  20 
. CONTINUE 

Evaluate  the  cubic. 

DXM=XX-X(I) 

DXP=X(I+1 )-XX 
DEL-X(I+1)-X(I) 

F=FDP( I ) ‘DXP. (DXP.DXP/DEL-DEL)/6 . 
k +FDP(I+1 )«DXM.(DXM.DXM/DEL-0EL)/6. 
I;  +Y(I).DXP/DEL+Y(I  + 1).DXM/DEL 

RETURN 

END 
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APPENDIX  F Equation  for  Combined  Seasonal  Efficiency  - EFFSPWT 


EFFSPWT  ; Efficiency  equotion 
Elo  • ox/(*+b)  + c 


September  22,  1988  Cheol  Pork 
Rev:  Septmeber  26,  1988 


PROGRAM  EFFSPWT 

REAL  EF,EFMULT. ETA225 , ETA225N , ETANEW , ETASS , ETAZERO , U | XDHW , 

k XLOAD,XSPWT 

REAL  ETA 

INTEGER  ] ,IMAX,J,JMAX.NCASE,NDATA 


PARAMETER  ( I MAX-50, JMAX-20) 

DIMENSION  XlOAD(IMAX) , ETANEW (I MAX. JMAX) .XDHW(JMAX) 
DIMENSION  U(JMAX) , EF( JMAX) . ETA225N( JMAX) 


OPEN(7. FILE-' EFFSPWT. IbP') 

OPEN( 8 .FILE-' EFFSPWT . OOT ' ) 

CLOSE(8. STATUS- 'DELETE') 

OPEN ( 8 . F I LE- ' EFFSPWT . OUT ' ) 

REWIND  7 

C Reod  the  input  dote  for  zero-voter- load. 

READ(7.*)  ETAZERO. ETA225. ETASS. EFMULT 

C Reod  the  number  of  cote*,  the  voter  drovn  per  day  in  gollon,  and 
C the  energy  foctor«(X).  ond  the  lood  foctor  due  to  domestic  hot  voter 
C drov(-). 

READ(7.*)  NCASE 
DO  10  J-1. NCASE 
READ(7.*)  U(J),E''(J).XDHW(J) 

EF(j)-EFMULT*EF(j)/iee 
ie  CONTINUE 

DO  20  J-1 .NCASE 
XSPWT-0.225+XDHW(J) 

ETA225N ( J )-ETA( ETAZERO . ETA225 . ETASS . XSPWT ) 

20  CONTINUE 

C Reod  the  loads  vhers  the  effciencies  ore  colculoted. 

REA0(7.*)  NDATA 
DO  30  1-1 .NDATA 
REAO(7.*)  XLOAD(l) 

XLOAD( I )-XLOAD( I )/100 
30  CONTINUE 

DO  40  J-1 .NCASE 
DO  40  1-1 .NDATA 

ETANEW( I . J )-ETA( EF( J ) . ETA225N( J ) . ETASS . XLOAD ( I ) ) 

40  CONTINUE 

DO  50  1-1 .NDATA 

WRITE(8.1000)  XLOAD(I).(ETANEW(l.J).J-1.NCASE) 

50  CONTINUE 

1000  FORMAT(0F10.4) 
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STOP 

END 


FUNCT I ON  ETA( ETAZERO . ETA225 . ETASS . XLOAD) 


REAL  A . B . C . ETA . ETA225 . ETASS . ETAZERO . X LOAD 

C A curve-fit  equotion 
OETAZERO 

B— 0 . 225*  ( ETASS-ETA225) 

k /(0  225* ETASS-ETA225+0. 775* ETAZERO) 

A-(ETASS-C)*(1 .+B) 

ET A=A . X LOAD/ ( X LO AD+B ) +C 

RETURN 

END 
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APPENDIX  G Efficiency  equation  and  bin  method  - EQTBIN 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 


EQTBIN:  Efficiency  equation  and  bin  analysis 

Eta  = ax/(x+b)  + c 


September  22,  1988  Cheol  Park 
Rev:  January  27,  1989 


PROGRAM  EQTBIN 


REAL  EF , ETA225 , ETA225N , ETANEW,  ETASS  , ETA2ER0 , U , 

& XDHW,XLOAD,XSPWT,ETAU,ETASBIN 

REAL  ETA 

INTEGER  I,IMAX,J,JMAX,NCASE,NDATA 
CHARACTER  TITLE*80 


PARAMETER  (IMAX=50 , JMAX=20) 

DIMENSION  XLOAD(IMAX),ETANEW(IMAX,JMAX),XDHW(JMAX) 

DIMENSION  U(JMAX)  ,EF(JMAX)  ,ETA225(JMAX)  , ETA225N(JMAX) 

DIMENSION  ETAU(IMAX),ETASBIN(JMAX) 

0PEN(7 , FILE=' EQTBIN. INP' ) 

OPEN ( 8 , FILE= ’ EQTBIN . ODT ' ) 

CLOS  E ( 8 , STATUS= ' DELETE ' ) 

OPEN ( 8 , FILE= ' EQTBIN . ODT ' ) 

REWIND  7 

C Read  the  input  data  for  zero-water- load. 

READ (7,*)  ETAZERO, ETASS 

C Read  the  number  of  cases,  the  water  drawn  per  day  in  gallon,  and 
C the  energy  factors(%),  and  the  load  factor  due  to  domestic  hot  water 
C draw(-). 

READ(7,*)  NCASE 
DO  10  J=l, NCASE 

READ(7,*)  U(J),EF(J),XDHW(J).ETA225(J) 

EF(J)=EF(J)/100. 

10  CONTINUE 


C Find  a new  combined  efficiency  at  the  point  where  the  space  load 

C is  22.5  percent. 

DO  20  J=l, NCASE 
XS PWT=0 . 2 2 5+XDHW ( J ) 
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ETA225N( J )=ETA(ETAZERO , ETA225 ( J ) , ETASS ,XSPWT) 

20  CONTINUE 

C Read  the  loads  where  the  effciencies  are  calculated. 

READ (7,*)  NDATA 
DO  30  1=1, NDATA 
READ(7,*)  XLOAD(I) 

XLOAD ( I ) =XLOAD ( I ) /I 00 
30  CONTINUE 

DO  50  J=1,NCASE 
DO  40  1=1, NDATA 

ETANEW(I,J)=ETA(EF(J) ,ETA225N(J) , ETASS, XLOAD(I)) 

IF(I  .GT.  1)  THEN 

ETAU(I-1)=  ETANEW(I,J) 

ENDIF 

40  CONTINUE 

CALL  ESEASONC  ETAU,  ETAS) 

ETASBIN(J)  =ETAS 

50  CONTINUE 

TITLE='  ******  OUTPUT  OF  EQTBIN.FTN  ******' 
WRITE(8,2000)  TITLE 
WRITE(8,3000)  (U(J) , J=1 .NCASE) 

WRITE (8, 4000) 

DO  60  1=1, NDATA 

WRITE(8 , 1000)  XLOAD(I) , (ETANEW(I , J) , J=1 , NCASE) 

60  CONTINUE 

WRITE(8 , 5000)  (ETASBIN(J) ,J=1 .NCASE) 

1000  FORMAT (12F1 0.4) 

2000  FORMAT (A80) 

3000  FORMATS  XLOAD  U=M2F10.1) 

4000  FORMAT  ( ) 

5000  FORMATC/'  Eta.bin  =M2F10.4) 


STOP 

END 


C 

C 

C 

C 

C 


FUNCTION  ETA(ETAZERO,ETA225, ETASS, XLOAD) 


REAL  A. B.C, ETA, ETA225, ETASS, ETAZERO, XLOAD 
C A curve- fit  equation 
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C=ETAZERO 

B=-0.225*(ETASS-ETA225) 

& /(O . 225*ETASS-ETA225+0 . 775*ETAZERO) 

A=(ETASS-C)*(1.+B) 

ETA=A*XLOAD/ (XLOAD+B) +C 

RETURN 

END 


C**icic**irk*irk*****icicirk***ic****irkic****icirk***irk**icirk*ic-k-k**-k*-kic-kic-k-k*****-*;*-kic**-k-k*ic-f; 

C 

c ESEASON  : seasonal  efficiency  for  combined  space  and  water  loads 

c 

c January  26,  1989  Cheol  Park 

c 


subroutine  eseason(etau,  etas) 

real  suml,  sum2,  xspace,  etau,  etas,  temp 

real  n 

integer  i,  imax 

parameter  ( imax=  15  ) 

dimension  xspace (imax) , etau(imax),  temp(imax), 

& n(imax) 

c bin  data  for  region  IV 

data  n/0.132,  0.111,  0.103,  0.093,  0.100,  0.109,  0.126, 

0.087,  0.055,  0.036,  0.026,  0.013,  0.006,  0.002, 

0.001/ 

data  temp/  62.,  57.,  52..  47.,  42.,  37.,  32.,  27.,  22.. 

& 17.,  12.,  7.,  2.,  -3.,  -8./ 

c design  condition 

data  tref/65.0/,  tdesign/5.0/,  oversize/0.7/ 


suml  =0.0 
sum2  = 0.0 

do  10  i=l,imax 

xspace(i)  = (tref-temp(i) )/( (tref-tdesign)*(l .+oversize) ) 

suml  = suml  + n(i)*xspace(i) 

sum2  = siim2  + n(i)*xspace(i)/etau(i) 

if  (i  .eq.  imax)  then 
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print  *,  'i,  xspace,  suml , suin2,  etau  =' , 1,  xspace(i)  .suml, 
& sum2,  etau(l) 

endif 

continue 

etas  = sunil/suin2 
print  * 

print  *,  ' Seasonal  Efficiency  = etas 

return 

end 
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